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Synopsis 



There is a growing consensus among physicists that the classical notion of spacetime has 
to be drastically revised in order to find a consistent formulation of quantum mechanics 
and gravity. One such nontrivial attempt comprises of replacing functions of continuous 
spacetime coordinates with functions over noncommutative algebra. Dynamics on such 
noncommutative spacetimes (noncommutative theories) are of great interest for a variety 
of reasons among the physicists. Additionally arguments combining quantum uncertain- 
ties with classical gravity provide an alternative motivation for their study, and it is 
hoped that these theories can provide a self-consistent deformation of ordinary quantum 
field theories at small distances, yielding non-locality, or create a framework for finite 
truncation of quantum field theories while preserving symmetries. 

In this thesis we study the gauge theories on noncommutative Moyal space. We find 
new static solitons and instantons in terms of the so-called generalized Bose operators 
(GBO). GBOs are constructed to describe reducible representation of the oscillator alge- 
bra. They create/annihilate k-quanta, k being a positive integer. We start with giving 
an alternative description to the already found static magnetic flux tube solutions of the 
noncommutative gauge theories in terms of GBOs. The Nielsen-Olesen vortex solutions 
found in terms of these operators also reduce to the ones known in the literature. On the 
other hand, we find a class of new instanton solutions which are unitarily inequivalent to 
the ones found from ADHM construction on noncommutative space. The charge of the 
instanton has a description in terms of the index representing the reducibility of the Fock 
space representation, i.e., k. After studying the static soliton solutions in noncommu- 
tative Minkowski space and the instanton solutions in noncommutative Euclidean space 
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we go on to study the implications of the time-space noncommutativity in Minkowski 
space. To understand it properly we study the time-dependent transitions of a forced 
harmonic oscillator in noncommutative 1+1 dimensional spacetime. We also provide an 
interpretation of our results in the context of non-linear quantum optics. We then shift 
to the so-called DSR theories which are related to a different kind of noncommutative 
(n-Minkowski) space. DSR (Doubly/Deformed Special Relativity) aims to search for an 
alternate relativistic theory which keeps a length/energy scale (the Planck scale) and a 
velocity scale (the speed of light scale) invariant. We study thermodynamics of an ideal 
gas in such a scenario. 

In first chapter we introduce the subjects of the noncommutative quantum theories 
and the DSR. Chapter 2 starts with describing the GBOs. They correspond to reducible 
representations of the harmonic oscillator algebra. We demonstrate their relevance in 
the construction of topologically non-trivial solutions in noncommutative gauge theories, 
focusing our attention to flux tubes, vortices, and instantons. Our method provides a 
simple new relation between the topological charge and the number of times the basic 
irreducible representation occurs in the reducible representation underlying the GEO. 
When used in conjunction with the noncommutative ADHM construction, we find that 
these new instantons are in general not unitarily equivalent to the ones currently known 
in literature. 

Chapter 3 studies the time dependent transitions of quantum forced harmonic oscilla- 
tor (QFHO) in noncommutative R 1 ' 1 perturbatively to linear order in the noncommuta- 
tivity 9. We show that the Poisson distribution gets modified, and that the vacuum state 
evolves into a "squeezed" state rather than a coherent state. The time evolutions of un- 
certainties in position and momentum in vacuum are also studied and imply interesting 
consequences for modelling nonlinear phenomena in quantum optics. 

In chapter 4 we study thermodynamics of an ideal gas in Doubly Special Relativity. 
We obtain a series solution for the partition function and derive thermodynamic quanti- 
ties. We observe that DSR thermodynamics is non-perturbative in the SR and massless 
limits. A stiffer equation of state is found. We conclude our results in the last chapter. 
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Notation and Abbreviations 



The metric for 1 + d Minkowski space has been taken to be 

r}^ = dia(l, -1, -1, ...) 

where the first element comes for time coordinate. 
The Abbreviations used in the thesis are 

• GBO: Generalized Bose Operator 

• UIR: Unitary Irreducible Representation 

• ADHM: Atiyah, Drinfeld, Hitchin and Manin 

• SD: Self Dual 

• ASD: Anti-self Dual 

• SHO: Simple Harmonic Oscillator 

• FHO: Forced Harmonic Oscillator 

• QFHO: Quantum Forced Harmonic Oscillator 

• SR: Special Relativity 

• DSR: Doubly /Deformed Special Relativity 
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Introduction 



There is a growing consensus among physicists that the classical notion of spacetime 
has to be drastically revised in order to find a consistent formulation combining quan- 
tum mechanics and gravity [1-3]. One such nontrivial attempt comprises of replacing 
functions of continuous spacetime coordinates with functions over noncommutative al- 
gebra [4, 5] . Dynamics on such noncommutative spacetime (noncommutative theories) 
are of interest for a variety of reasons among the physicists for a long time. Let us start 
with introducing such noncommutative theories. 

1.1 Noncommutative theory 

The idea of extension of noncommutativity to the coordinates was first suggested by 
Heisenberg as a possible solution for removing the infinite quantities of field theories 
before the renormalization procedure was developed and had gained acceptance. They 
seem to have first appeared in a letter from Heisenberg to Peierls in 1930 [6]. The first 
paper on the subject was published in 1947 by Hartland Snyder who sought to use it to 
regularize Quantum Field Theories [7]. The success of renormalization theory however 
drained interest from the subject for some time. 

The renewed interest by particle physics community started to grow after the paper 
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by Seiberg and Witten [8]. Arguments combining quantum uncertainties with classical 
gravity also provide an alternative motivation for their study [1], and these theories 
can provide a self-consistent deformation of ordinary quantum field theories at small 
distances, yielding non-locality [9-11], or create a framework for finite truncation of 
quantum field theories while preserving symmetries [12-17]. Also, the emergence of 
noncommutative field theories in string theory [8, 18, 19] has provided a considerable 
impetus to their investigation. 

The noncommutativity of the coordinates can be described by the following commu- 
tation relation 



where 6^ v is an anti-symmetric tensor. The simplest case corresponds to 6^ u being 
constant. The field theories on such spacetimes are described by elements of unital 
algebra (associative, but not commutative) generated by x^s. There is a sense in which 
the algebra over the noncommutative space R# and the algebra over commutative space 
Ro (d being the dimension of the space) have the same topology. This algebra can 
be considered as a deformation of the algebra over Rq such that the elements of the 
algebra have the same addition law, but a different (noncommutative) multiplication 
law which reduces to the commutative point-wise multiplication in the limit 6^ — > 0. 
This noncommutative multiplication law is often denoted by / * g or the star product 
to distinguish it from the ordinary pointwise multiplication of functions. Let us give an 
example of such star product. If we denote by m the pointwise multiplication map 



x u ~\ i0 l 



(1.1) 




(1.2) 



then the star product is generally given by 



(/ * g) (x) =m e (f <g> g) (x) =m oj 9 (/® g) (x) 



(1.3) 
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For instance, Fq corresponding to the Moyal product is given by 

J e — e k e ^{ dx ^® dx " . (1.4) 

One can easily check that the following is satisfied: 
1.1.1 The Poincare invariance 

Noncommutativity of the spacetime coordinates apparently conflicts with Poincare in- 
variance, as one can see in (1.1). We circumvent this problem by "twisting" the Poincare 
symmetry [20]. It goes as follows. First of all one needs to define the action of a Poincare 
group element on a tensor product space, which is generally given by the coproduct. The 
usual coproduct used is as follows: 

A (A)=A®A, (1.6) 

A being an element of the Poincare group. We twist the coproduct to 

Ae(A) = ^ 1 A (A) T 6 = F e \A ® K)T e (1.7) 

Action of a Poincare element on star product of two functions is given by 

A> (f * g) = A> m e (f ® g) = m e A e (A) (f ® g) (1.8) 

Putting f,g = x M , x v one can check that the commutation relation (1.5) remains invariant 
under the action of Poincare group. This modification in coproduct changes the standard 
Hopf algebra structure associated with the Poincare group (the Poincare-Hopf algebra) 
to the twisted Poincare-Hopf algebra. This twist in turn give rise to twisted statistics 
leading to variety of interesting consequences [21-31]. We will not be discussing them 
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here to avoid any deviation from the desired topic. 

In this thesis we study the gauge theories on noncommutative Moyal space 
{Qjiv = constant). We start by studying the implications of the construction of so-called 
generalized Bose operators (GBO) [32] in noncommutative gauge theories. GBOs corre- 
spond to the reducible representation of the oscillator algebra 



We discuss them in detail in chapter 2. The noncommutative theory discussed here 
contains commutative time. Having studied gauge theories with space-space noncommu- 
tativity through the construction of GBOs we proceed to understand the implications 
of noncommutative time coordinate. There have been even claims that quantum field 
theories based on noncommutative time coordinate are nonunitary [33]. In contrast, in a 
series of fundamental papers, Doplicher et al. [1,34] have studied (1.1) in complete gener- 
ality, without assuming that 9 0i = and developed unitary quantum field theories which 
are ultraviolet finite to all orders. In this thesis we study a simple quantum mechanical 
model of forced harmonic oscillator in 1 + 1 dimensional noncommutative spacetime to 
get the essence of the effect of noncommutative time coordinate. For this we follow the 
unitary formulation of quantum mechanics with time-space noncommutativity discussed 



There are kinds of noncommutative spaces other than the Moyal space. In particu- 
lar, we study the thermodynamics in the so-called Doubly/Deformed special relativistic 
(DSR) theories which are related to a noncommutative space known as /t-Minkowski 
space [36-41]. DSR (Doubly/Deformed Special Relativity) aims to search for an al- 
ternate relativistic theory which keeps a length/energy scale (the Planck scale) and a 
velocity scale (the speed of light scale) invariant. 



[a, a)] = 1. 



(1.9) 



in [35]. 
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1.2 Doubly /Deformed Special Relativity 

Attempts to combine gravity with quantum mechanics in search of the theory for quan- 
tum gravity always seem to give rise to the Planck length (l P = \f^j that provides 
the scale at which the quantum effects of gravity will show up [42-52]. The existence 
of such a length scale is in conflict with the equivalence principle because observers in 
different inertial frames will not agree on lp due to the Lorentz-Fitzgerald contraction. 
To understand in a better way consider the following: If the Planck scale is to play the 
role of the threshold for the discrete spacetime, it may come in the quantities involving 
the metric (for example the action). For equivalence principle to hold, i.e., for physics 
to be independent of the inertial frames the action and hence the Planck scale has to 
remain invariant under the relativistic transformation. 

It has been shown that it is possible to still have equivalence principle by deforming 
Special Relativity (SR). These classes of theories fall under the name Doubly/Deformed 
Special Relativity (DSR) [53-56]. In DSR, apart from the constancy of speed-of-light 
scale, the Planck length lp or equivalently Planck energy k is also constant under coor- 
dinate transformation from one inertial frame to another. It was first proposed by G. 
Amelino-Camelia [53] . To proceed one can take lessons from the transition from Galilean 
to Lorentz transformations. The Galilean transformation is linear in velocity and hence 
it does not keep any velocity scale invariant. To have an invariant velocity scale we 
needed to accept the Lorentz transformation which is nonlinear in velocity. Thus if one 
wants to have Planck length/energy scale invariant the relativistic transformation has 
to be nonlinear in position/momentum space. One also keep in mind the possibility of 
the symmetry group being different from the usual Lorentz group. This symmetry group 
must act nonlinearly in the position/momentum space representation to keep Planck 
length/energy invariant. In the special relativistic limit, Planck length/energy — > O/oo 
the symmetry group has to go to the Lorentz group. 

J. Magueijo and L. Smolin have argued that the symmetry group for the DSR trans- 
formations remains to be the Lorentz group [57,58]. But in DSR prescription they act 
nonlinearly on the position/momentum space. They have also suggested the following 



Chapter 1. Introduction 



6 



l + (7"l)f- 


7 u ^f 


l(Pz - ve] 


1 


l + (7-l)f" 


7 v ^f 


Px 




! + (7-l)f- 


7 v ^f 


Py 




l + (7-l)^- 


7^ 



nonlinear representation of the Lorentz transformation on momentum space 

a(-)-^' = ; { i(£ ~: Pz) P (i-io) 

A(«):p.^Pi - 1 + ( ; ( -^- 7 ^ a- 11 ) 
A(t,) : p x j/ x = - ; P ;, £ (1.12) 

(1.13) 

with 7 = ^7=f,c = 1. The above transformation corresponds to the boost along z-axis. 
Note that they also satisfy the group multiplication law of the usual Lorentz group, i.e., 

M«)A(«)=A(££) (1.14) 
The transformation keeps the following mass invariant 

m= l r^— 2 (1-15) 

We call it the invariant mass. Note that unlike the special relativistic (SR) case, this is 
not same as the rest mass energy. 

There have been attempts to study the representation on the position space too 
[59-61]. But there is no proper understanding of its relation to the momentum space 
representation. In this thesis we have used the momentum space representation given 
by (1.10) - (1.13) which in turn implies the modification in the dispersion relation 

£ 2 -p 2 = m 2 - ^j 2 . (1.16) 

Consequences of the modified dispersion relations on the thermodynamics are being 
studied extensively to infer the effect of Planck scale physics [62-68]. The effect of 
modified dispersion relations in loop-quantum-gravity on black hole thermodynamics was 
studied in [63]. The same as a Lorentz violating phenomena on the thermodynamics of 
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macroscopic systems (like white dwarfs) [64,66,67] and as a noncommutative phenomena 
on cosmology and astrophysical systems [62,68] have also been studied. Moreover, photon 
gas thermodynamics in the context of modified dispersion relations [65] and DSR [69] are 
being investigated. In [69] the effect comes solely because of the presence of a maximum 
energy scale as the photon dispersion relation remains unmodified. In this thesis we 
study the thermodynamics of an ideal gas in the above described DSR scenario. 

The organization of the thesis is as follows: Chapter 2 starts with describing the 
GBOs. They correspond to reducible representations of the harmonic oscillator algebra. 
We demonstrate their relevance in the construction of topologically non-trivial solutions 
in noncommutative gauge theories, focusing our attention to flux tubes, vortices, and 
instantons. Our method provides a simple new relation between the topological charge 
and the number of times the basic irreducible representation occurs in the reducible 
representation underlying the GBO. When used in conjunction with the noncommutative 
ADHM construction, we find that these new instantons are in general not unitarily 
equivalent to the ones currently known in literature. 

Chapter 3 studies the time dependent transitions of quantum forced harmonic oscilla- 
tor (QFHO) in noncommutative R 1 ' 1 perturbatively to linear order in the noncommuta- 
tivity 9. We show that the Poisson distribution gets modified, and that the vacuum state 
evolves into a "squeezed" state rather than a coherent state. The time evolutions of un- 
certainties in position and momentum in vacuum are also studied and imply interesting 
consequences for modelling nonlinear phenomena in quantum optics. 

In chapter 4 we study thermodynamics of an ideal gas in Doubly Special Relativ- 
ity. New type of special functions (which we call Incomplete Modified Bessel functions) 
emerge. We obtain a series solution for the partition function and derive thermodynamic 
quantities. We observe that DSR thermodynamics is non-perturbative in the SR and 
massless limits. The equation of state found makes the pressure vs energy density graph 
stiffer. 

We conclude our results in the last chapter. We take this opportunity to mention 
that we have corrected some typos in the notations and some results in the published 
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version of our papers [70] and [71]. While the other corrections are trivial (or they come 
in the middle of calculations) we must mention the corrected results in (3.140)-(3.142). 



Chapter 2 

Noncommutative Vortices and 
Instantons from GBOs 

2.1 Introduction 

Detailed investigations of noncommutative gauge theories have led to the discovery of 
localized static classical solutions in noncommutative spaces [72-80]. Among the models 
of gauge theories in noncommutative spaces, one of the simplest is the Abelian-Higgs 
model possessing vortex-like solutions [81-84]. Another interesting class of solutions in 
noncommutative euclidean space are instantons in U (N) Yang-Mills theories. Nekrasov 
and Schwarz developed a generalization of ADHM construction as given in [85] to find 
these noncommutative instantons [86]. The U(N) Yang-Mills instantons in R^c* an d 
Ryvc x R^, were further studied in [87-93]. Pedagogical reviews can be found in [94,95]. 
Apart from these, there are multitudes of other solutions in noncommutative gauge fields 
like merons [93], flux tubes [96], monopoles [97], dyons [98], skyrmions [99], false vacuum 
bubbles [82], to name just a few. 

In this chapter we present a new construction of such topological objects, based on 
an analysis of the reducible representations of the standard harmonic oscillator algebra. 
Our method gives rise to new instanton solutions (i.e. not gauge equivalent to the known 
ADHM instantons) on a noncommutative space, and in the process provides a simple 
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interpretation for the instanton number: it simply "counts" the number of copies of the 
basic irreducible representation. 

Our construction relies on operators called Generalized Bose Operators (GBO) [32, 
100] which provide an explicit realization of the reducible representations of the oscillator 
algebra, and are well-known in the quantum optics literature (see for examples [101,102]). 
As a warm-up, we first study the significance of GBOs in constructing fluxes and vortices 
with higher winding numbers, and then discuss instanton solutions in noncommutative 
Yang-Mills theories. 

The chapter is organized as follows. We start with a brief review of GBOs and their 
representations in section 2.2. In section 2.3, we discuss the flux tube solutions of [96] 
in the language of GBOs and then we go on to show the relevance of these operators 
in noncommutative Nielsen-Olesen vortices. Section 2.4 discusses the noncommutative 
instantons. Using the GBO in conjunction with the ADHM construction, we construct a 
class of new instantons and compute their topological charges. Notations and formalism 
of the theories discussed in sections 2.3.1, 2.3.2 and 2.4 are all independent to each other 
and they should not be mixed. 

2.2 Generalized Bose Operators - A Brief Review 

Brandt and Greenberg [32] give a construction of Generalized Bose Operators that change 
the number of quanta of the standard Bose operator a by 2 (or more generally by a 
positive integer k). We briefly recall their construction in this section. 

Consider the infinite-dimensional Hilbert space T-L spanned by a complete orthonormal 
basis {\n), n — 0, 1, • • • , oo} labeled by a non-negative integer n. Vectors in TL are of the 
form = J2 n c n\ n )i c n £ C Vn such that J2 n \ c n\ 2 < oo. We formally write 




(2.1) 
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The standard bosonic annihilation operator a acts on this basis as 

a\n) = n*\n- 1), Vn > 1 and a|0) = 0. (2.2) 
The annihilation operator is unbounded, and hence comes with a domain of definition: 

n|c n | 2 < oo \ . (2.3) 



{ n n 



Its adjoint a* satisfies 

a t|n) = (n + l)s|n + 1), Vn > (2.4) 
and has the same domain D a . The number operator 

N = at a (2.5) 

is defined with a domain 



Vn = } ^2 c ™\ n ) : ^2 n2 \ c n\ 2 < oo > . 

I n n ) 



(2.6) 



The basis vectors {\n}} are eigenstates of N: 

N\n)=n\n). (2.7) 
On Vn, the operators a and a* satisfy 

[a,a f ] = l. (2.8) 

While iV counts the number of quanta in a state, the a and a* destroy and create 
respectively a single quantum. They satisfy the following commutation relations 

[N,a] = -a, [N,<J]=a! (2.9) 
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We denote this irreducible representation of the oscillator algebra by (a, %). This, up to 
a unitary equivalence, is the unique irreducible representation [103]. The Hilbert space 
% can be split into two disjoint subspaces % + and T-L- given by 

H + = {X>2n|2n> EH) (2.10) 
U - = {$> 2n+1 |2n+l> G^} (2-11) 

such that 

H = H + ®H^ (2.12) 

The projection operators 

A+ = |2n)(2n| = cos 2 (^f), and A_ = ^ \2n + l)(2n + 1| = sin 2 (^~) 

n=0 ^ ' n=0 ^ ' 

(2.13) 

project onto the subspaces T-L + and "H_ respectively. One can define the operators b± 
and its adjoint b± on the subspaces %± by 

b + \2n) =n5|2n-2), ^|2n) = (n + 1)5 \2n + 2), 6+|0) = (2.14) 

6_|2n + l) =n5|2n-l), &L|2n + 1) = (n + l)^\2n + 3), 6_|1) = (2.15) 
with domain of closure V a fl H±. On the domain T>n fl "H± we have 

[k,&L] = l- (2-16) 

Thus (&_,%_), (b +: H + ) and (a, "H) are isomorphic to each other. In other words, there 
exist unitary operators U± such that 

[/ ± 6 ± [4 = a. (2.17) 
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Using the projection operators A±, one can define an operator b as 

b = b + A+ + b_A_ (2.18) 

on H whose action on the basis vectors \n) is 

b\2n) =ns|2n-2), b\2n + 1) = rai \2n - 1). (2.19) 

Notice that both |0) and |1) are annihilated by b. The operator b satisfies the commuta- 
tion relation 

[TV, 6] = -26. (2.20) 

The adjoint of b satisfies 

tf\2n) = (n + l)s|2ra-2), tf\2n + 1) = (n + 1)3 \2n - 1). (2.21) 

and 

[N,b jf ]=2tf. (2.22) 
A new number operator can be defined as 

M = b ] b= -(N-A-) (2.23) 

which has the same eigenstates |n), but each eigenvalue is two-fold degenerate! We 
denote these eigenvalues by 

m n = ^(n-A n _) (2.24) 

where we define \ n ± as 

1 if n is even 

(2.25) 

if n is odd 



A 



n+ 



(n|A + |n) 



2 ( nn \ 

cos It 
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n /mr\ if n is even 
A n _ = (n|A_|n) = sin 2 (— ) = { (2.26) 

v 2 y [ 1 if n is odd 

We rewrite the action of b and in a compact form 

b\n) = ml\n - 2) and &t|n) = ( m n + l)*\n + 2). (2.27) 

The operators b and have £> a as their domain of closure and they satisfy 

[6,61 = 1 (2.28) 

in the domain T>n- Hence (b, ft) forms a reducible representation of the oscillator algebra 
characterized by (2.8) having (2.18) as its irreducible decomposition. 

One can generalize the above construction and formulate an operator which lowers 
a state |n) by k— steps. First of all we note that ft can be split as 

ft = Ut ] © ft? © ... © Ut\ (2.29) 

for an integer k where ft^'s are defined as 

Uf = \^c kn+l \kn + i) e ft j . (2.30) 
We define projection operators A\® as 

oo 

Af } = ^2\kn + i)(kn + i\, % = 0, 1, • • • k - 1. (2.31) 

n=0 

which project onto the subspace ft^. In each subspace ftf^, one can define operators 
b\® and their adjoints b\ k ^ that satisfy 

bf\bf )] ] =1 (2.32) 
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and hence correspond to the UIR of the oscillator algebra. They act on the states as 



bf^\kn + i) — y/n\kn + i — k), b\ k ^\kn + i) — \Jn + l\kn + i + k) 



(2.33) 



We can easily check that 



6? ) >A (*) = 6 wt A w 



A reducible representation is given by 



k-l 



k-l 



(k) 



i=0 



They satisfy the oscillator algebra commutation relation 



(2.34) 



(2.35) 



[&(*),&(*)t] = i 



(2.36) 



on the domain T>^. Thus {bf^\l-if^) for i — 0, 1, k — 1 are isomorphic to (a,T-L) and 
hence (6^ fc - ) ,'H) forms a reducible representation of the oscillator algebra. Discussions 
from (2.10) to (2.28) represent the case k = 2, the simplest non-trivial example of this 
construction. Henceforth we will use b for b^. An explicit expression for b is [100] 



1 ( 1 A 1 

o = I a— aA + + a 



\/2 V VA 7 



(2.37) 



Before we end this section, let us point out a minor generalization of the Brandt- 
Greenberg construction. Under any unitary transformation U± defined on %± that 
transforms b± as b± — > U±b±U±, the fundamental commutation relation (2.16) remains 
unchanged. In particular, if we choose 



U ± (z ± )=e z ^ 



z±b\-z±b± 



(2.38) 



then we find that 



b ± (z ± ) = U±(z±)b±Ul(z±) = 6± - z ± , 



(2.39) 
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i.e., we get the "translated" annihilation operator. One can construct a reducible repre- 
sentation using b + (z+) and b-(z-) as 

b(z+, z-) = b + {z + )A + + b^{z„)A_ = 6+A+ + 6_A_ - z + A + - z_A_. (2.40) 

Here b gets translated by different amounts in different subspaces H± and the "trans- 
lated" operator b(z + ,Z-) is unitarily related to b as 

b(z+, z.) = U(z+, z-)btf(z+, z-), U(z+, z.) = U+(z+)A + + )A_. (2.41) 



More generally, using (2.35) we can write 



k-l 

b^-J2 z M k) (2-42) 

i=0 



fc-i 

Y.UM)^- (2-43) 

i=0 



& (fc) (*o»2i» • • • ,z fc _i) = 
and the unitary operator is 

U (zq, Zi, ■ • • , Zk-\) = 



Though minor, this generalization will play a role in the construction of noncommutative 
multi-instantons. 

There exist other possibilities as well. For example, choosing 

£4(fc) = e^-#> (2.44) 

gives 

&±(£±) = 6± cosh^i + bl sinh£ ± . (2.45) 



The above is the well known squeezed annihilation operator. A reducible representation 
may be constructed: 

= M£+)A+ + M£-)A- (2.46) 
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and more generally, 

fe-i 

6(&,ei,---,&-i) = E 6 i fc) te) A 5 fc) ( 2 - 47 ) 

We discussed this "squeezing" operators just to give an example and will not study it in 
detail as they are not directly used further. 

Fields in noncommutative Moyal space are generally interpreted in terms of the simple 
harmonic oscillator algebra elements. Having found the reducible representation of this 
algebra in terms of the so-called GBOs, we try to seek nontrivial solutions of field theories 
on noncommutative spaces in terms of these operators. 

2.3 Static Solutions In Noncommutative Gauge The- 
ories 

As we are interested in exploring the relevance of the GBOs in noncommutative field 
theories, we start with the following two simple situations where this relevance is most 
visible: 

• the flux tube solution in (3 + 1)— dimensional pure gauge theory [96] 

• the vortex solution in (2 + 1)— dimensional abelian Higgs model [81]. 

We find a new interpretation for the already known solutions of the above theories in 
terms of the GBOs. 

2.3.1 Flux Tube Solution In Noncommutative Gauge Theories 

Consider pure U(l) gauge theory in (3 + l)-dimensional spacetime with only spatial 
noncommutativity. This theory incorporates magnetic flux tube solutions [96] which are 
important in the context of monopoles and strings discussed in [97]. We will search for 
non-trivial solutions of the static equation of motion. These solutions do not possess a 
smooth 9 — > limit [96], implying that they have no commutative counterpart, i.e, the 
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origin of this effect is entirely due the noncommutativity of the underlying space. The 
noncommutativity is only among the space coordiantes (time is commutative) 

[x\x j ] = i9 ij ; i,j = 1,2,3. (2.48) 

The anti-symmetry of the real ^-matrix, i.e., 

e ij = -e ji (2.49) 

guarantees a choice of axes in which the noncommutativity becomes 

[x 1 , x 2 } = i6, [x\x 3 ] = 0, [x 2 , x 3 } = (2.50) 

so that only the x 1 — x 2 plane is noncommutative. On a noncommutative space, "func- 
tions" are elements of the noncommutative algebra generated by the operators x l . Deriva- 
tives in the x 1 and x 2 directions are defined via the adjoint action 

d xi f='-[x 2 ,f], d x ,f = --[x\f] (2.51) 

while the derivatives in the x 3 and t directions are the same as in the commutative 
case. We can define a set of complex (noncommuting) variables z and z and a set of 
creation-annihilation operators as 

z = ^(x 1 +ix 2 ), z = -^(x 1 -ix 2 ), a=-^=z, a} = -^=z (2.52) 

Here a and a* satisfy (2.8). With this convention, the derivatives with respect to the 
complex coordinates are given as 

d M f = -±=[a*,f], d s f=^=[a,f\. (2.53) 
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Integration on x 1 — x 2 plane is replaced by trace over Fock space H: 

/oo 
dx 1 dx 2 f(x\x 2 ) ->■ 2n9Ti n f(x\x 2 ) = 2n9 ^(nlftx 1 , x 2 )\n). (2.54) 



n=0 



where |n)'s are the number eigenstates and the factor 2n9 ensures the proper commuta- 
tive limit (6->0). 

We will follow the construction of gauge theories on a noncommutative space as given 
in [104]. In accordance with the commutative theory the field strength is defined as 

F lw = [V li ,V v \ (2.55) 

where the noncommutative covariant derivative is defined as 

V = d t + A , V 3 = d x3 +A 3 , V 1 = l -£ 2 + A 1 , V 2 = - l -x l + A 2 . (2.56) 

u u 

Here A^ is an anti-hermitian operator field 

K = -A, (2.57) 

The gauge transformation is such that it transforms £> M covariantly 

T? M ->■ UV^U' 1 (2.58) 

Noncommutative gauge theory has an alternative (and possibly more natural) formula- 
tion in terms of rather than A^. In terms of V^, for a pure gauge field, one works 
with the action 



S = ~ I dx'dtTriF 2 ,} (2.59) 
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This formulation of noncommutative gauge theory is classically equivalent to the stan- 
dard formulation where action is given by 



7l6 f 



dx dt Tr{F^ v } (2.60) 

y j 

with standard field strength given by 

= d,A v - d u A, + [Ap, A v \. (2.61) 
The derivatives for //, v — 1,2 are defined in (2.51). One can verify 

F^v = F^v + ^V, — -Q e otivz- (2.62) 
The anti-symmetric Levi-civita tensor is defined as 

1 if fiupcr is in even permutation of 0123 

— 1 if fxupa is in odd permutation of 0123 ; e^po- = —e iivpc '. (2.63) 
otherwise 

This also implies 

uv = -uT (2.64) 
Now it is just a matter of simple algebra to see that 

S = S+^- [ dtdx 3 Tr (1 + 2idF l2 ) . (2.65) 
9 J 

While each term of F 12 is a commutator or a total derivative and hence will contribute 

only to boundary terms, the term with Tr 1 is only an infinite irrelevant term. As the 
variations in X> M and A^ are same, both S and S give the same equations of motion and 
hence are classically equivalent! 

For static, magnetic configurations (d t = = A ) and with the choice d^Ai = 0, A 3 = 
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0, the equations of motion in terms of Ws reduce to 

[V,[V,V]]=0, with V= 4=(Z>i + iZ> 2 ), V=^ r (V 1 -iV 2 ) = -Vy (2.66) 

v 2 v 2 

It is easy to check that the standard "vacuum" configuration corresponding to A\ = 
A 2 = (and hence = 0) will satisfy the equation of motion (2.66) and corresponds 
to 

V = W (2 ' 67) 

We can construct solutions about this vacuum by taking a rotationally invariant ansatz 
V = af(N) and it can be shown that there exists a solution of the form (see Appendix 
A. 1 for a detailed calculation) 



P =^V^V^Xi |ri) H N = a\ n = 0,1, 2.... (2.68) 



n=no 



with N being the number operator as given in (2.5). This solution corresponds to 

oo 
n=rio 

The standard field strength becomes 

= ^ - = -u^v E l n )( n l- ( 2 - 70 ) 

?1=0 

Thus it represents a classical localized static circular magnetic flux tube in x 3 direction 
centred about origin of the (x 1 , x 2 ) plane with no related to its radius. The total magnetic 
flux $ gets quantized 

$ = 2ir9Tr(iF 12 ) = -2nn . (2.71) 
The choice n = corresponds to the vacuum configuration and has zero magnetic flux. 
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By virtue of the relation (2.28), 



V=JL (, 72) 



is also a vacuum solution (though reducible) of (2.66). Again, we can start with the 
ansatz V = bG(N) to construct the following solution (see Appendix A.l for a detailed 
calculation) 



v = ^- »^- »o-A- , £ |n)( „ |A+ + £ |n)(n|A _ , , ( , 73) 



\n=no+ n=no- 



Here n + = 0, 2, 4, ... and n _ = 1, 3, 5, .... This solution can be re-written in the form 

oo 

V = ^2 g{n)\n)(n + 2| (2.74) 



n=0 



with 



#(2n) = x /n - n' 0+ + 1 for n > n' 0+ , n' 0+ = ^± (2.75) 



#(2n + 1) = - n' _ + 1 for n > n' _, n „ = n ° - . (2.76) 

This solution is same as the higher moment solution obtained in [96] starting with the 
ansatz V = F(N)a 2 . Again, for the choice n + = and n _ = 1, this solution reduces 
to (2.72). 

Furthermore, the rotationally invariant ansatz of V — f(n)b^ about the reducible 
"vacuum" 

b (k) 

V^^ (2.77) 
for k > 2 gives the following solution of the equation of motion 



fe-i 

* k— 1 oo 



v 2=0 n=noi 
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with n 0i = kn + i, n — 0, 1, 2, .... These solutions are also same as the higher moment 
solutions obtained from the ansatz V = g(n)a k in [96]. Further computation of the 
magnetic field and other things for the above solution can be found in [96]. These 
solutions with the GBO represent static magnetic flux tubes with localized flux, with 
size of the tube in the i-th irreducible part of the Fock space related to n 0i . Thus the 
radial profile of the magnetic field is determined by the set {noi}. 

Let us try to understand the above construction in the following manner: The flux 
tube solution using irreducible representation (a, H) can be described by a single integral 
index n which is related to the extension of the magnetic field in the Fock space (see 
(2.70)). The use of reducible representation (b ( ^ k \'H) separates the Fock space in /c-parts, 
resulting in k numbers of indices n 0i each of which is related to the extension of magnetic 
field in the corresponding subspace of the Fock space. Increase in the number of indices 
(charges) clearly expand the moduli space of the static magnetic flux tube solutions! 
These solutions are nonperturbative in 9 as the reducibility of the Fock space and hence 
the notion of GBOs are typically attributed to noncommutative spaces! 

In this section we saw that there already exist certain solutions of noncommutative 
gauge theory which can be re-written in terms of the GBOs. This fact shows the impor- 
tance of GBOs and motivates us to seek such solutions in other noncommutative gauge 
theories. 

2.3.2 Noncommuative Abelian Higgs Model - Nielsen-Olesen 
Vortex Solution 

The abelian Higgs model in noncommutative spaces is of some interest because of its sim- 
plicity as a noncommutative gauge theory and the existence of vortex solutions. Various 
topologically non-trivial vortex solutions in this context have been studied in detail. 

An interesting class of vortex solutions in this theory is studied in [81], which are 
analogous to the Nielsen-Olesen Vortices in the commutative space [105]. The model 
is in (2 + 1)— dimensions, and consists of a complex Higgs field <3> which is a left gauge 
module (the gauge fields multiply the complex Higgs field <£> from left and <£> from right). 
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Minimizing the static noncommutative energy functional, Bogomolnyi equations are gen- 
eralized to the noncommutative space and 1 /^-expansion is done in large 6 limit. The 
equations are then solved order by order and the corrections to the leading order equa- 
tion converge rapidly. In the large distance limit (which is the commutative limit in 
this case), the solution reduces to the Nielsen-Olesen vortex solution in ordinary (com- 
mutative) abelian Higgs model. Here we will give a brief formalism of the standard 
noncommutative theory. For extensive calculation one can refer to [81]. 

The noncommutativity is same as that in section (2.3.1), with the only difference 
that now the space is 2— dimensional, so the direction x 3 is absent: 

[x\x 2 ]=i0. (2.79) 

The notations defined in (2.51), (2.52), (2.53) and (2.54) remain intact. The energy 
functional in the static configuration is given by [81] 



£ = Tr 



-(B + - l) 2 + V 2 <5>V Z <5> + T 
2 



(2.80) 



where V is a covariant derivative with the gauge field A and is given by 



V z $ = d z $ - iA z $, V z § = d z $ - iA z § 

V z § = d z § + i$A z , V z § = d-3 + i$A- z . (2.81) 

Note the left and right actions of the gauge fields on <3> and 4> respectively. The magnetic 
field B is defined as 

B = -t(d z A- z - d- z A z ) - [A z , A- z \ (2.82) 

With the abuse of notation a dimensionful quantity (3 — ^ (e and A are the coefficients 
of gauge coupling and the self-coupling of the Higgs fields respectively) has been set to 
1. For more details see [81]. T is the topological term defined as 

T = d m S m -i[A m ,S m ] + B. (2.83) 
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where m takes the values z and z. Here 

S m = -e mn ($£> n <J - £>„$$). (2.84) 

with the convention e zz = 1. It can be shown that TrT corresponds to the topological 
charge. Our prime interest is to study the Bogomolnyi equations. Minimizing (2.80) we 
get the following operator equations: 

V- z <& = 0, V z § = 0, £ = 1 - $$. (2.85) 

which are the noncommutative Bogomolnyi equations. Now one can do a 1/9 expansion 
of the Higgs and the gauge fields, in the large 9 limit 

$ = $oo + + 

9 (2.86) 

= ^(Ax> + g^-l + ) 

^ s = V$( 00 l^- 1 ) 

The factor of ^ is used for scaling the variables A as it is a 1-form. This makes sure 
that the fields $oo and are dimensionless. From the definition of the magnetic field 
we see 

£ = ^(£oo + ^£-i + ) (2.87) 

with 

= i([a,A 0O ]+ [a^An]) - [4x>,4x>] (2.88) 
£-1 = i([a,A_ 1 ]+[a t ,A_ 1 ])-([A 0O ,A_ 1 ] + [A_ 1 ,A 0O ]). (2.89) 

With this expansion, we can get the leading order 0(0) Bogomolnyi equation as 



(2.90) 
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This equation admits a solution [106] 

$oc = ~^=^a n , = ^ n -=L= (2.91) 

which represents an n-vortex at origin. A more general solution is discussed in [81] which 
represents n single vortices at n different points in the noncommutative plane and (2.91) 
is a special case of that general solution. But for our discussion, (2.91) is sufficient and 
due to its simple form, computation and understanding becomes easier. For the solution 
(2.91) we can check 

oo 

$oo$oo = ^|™>H^l- (2-92) 

m=n 

The next order Bogomolnyi equations become 

[a,$oo] = ^oo$oo, [a\$oc] = i&ocAoc (2.93) 
which can be solved to get (for details see appendix A. 2) 

= -i - 1 a (VN - VN + n) , A 00 = i (Vn - VN + n) a f - 1 (2.94) 
VATT V / V ) y/N + 1 

N being the number operator (2.5). 

In the coherent state \u>) (a\u>) = ou\u>)), the expectation of the field $oo is 

M$oola;) = u n (uj\ 1 \u) with u = \oo\e i<p . (2.95) 
\Ja n a) n 

The phase dependence is e m,p , which comes solely from u n as the other factor (u \ ^=f^ 
\u) is purely real, signifying a vortex in the noncommutative plane. The large distance 
behavior is given by the large u limit or equivalently large (N) limit [81]. The coherent 
state expectations in this limit becomes (for details see appendix A. 3) 

— n n 

(wl^oolw) « e mip , (wlAoolw) « i — , (w|4x>|w) « -i — , (2.96) 
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which is exactly like the commutative n-Nielsen-Olesen Vortex. 

It is interesting to note that the leading order magnetic field is (for details see A. 4) 

B oo =n\0)(0\ (2.97) 

which means that the magnetic field of the solution is localized and the magnetic fluxes 
are confined. The flux (trace of the field) is quantized and is characterized by the integer 
n. 

We try to seek the vortex solutions in terms of the GBOs. One can easily check that 

^new = 1 h n $new = &t n }_ ( 2 .98) 

satisfies (2.90). This gives (details in appendix A. 2) 

A new = _i( a _ _^ =b n ab \n * A ™ = V f _ _^ =& n fl t & tn * ^ _ 

V. Vb n tf n Vb n tf n J V vW" Vb n tf n J 

(2.99) 

The new solutions satisfy 

oo 

$oo$oo = ^ |m)(m| + I. (2.100) 

m=2n 

The expectation value in the coherent state \oj) (eigenstate of a) gives a phase dependence 
of e l2n<p (as can always be reduced to the form F(N)a 2n ), a characteristic feature 
of 2n vortex in noncommutative plane. In the large u limit it gives the large distance 
behavior: 

— 77 77 

{u\$™\u) « e i2nv , « i-, (u\A™ w \uj) « -i- (2.101) 

which is exactly the commutative 2n Nielsen-Olesen vortex. One can also calculate the 
magnetic field for the new solution to be (see Appendix A. 4). 

B™ w = 2n|0)(0|. (2.102) 
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Till now all the calculated properties of the new solutions matched exactly with those 
of the Witten's solution (2.91) of vortex number In. This stimulates us to compare the 
expression of the new solution (2.98) with that of the Witten's solution. For the simplicity 
of expression and better understanding of the underlying algebra, we take n — 1 in (2.98). 
Using the explicit expression of the GBO b, the new vortex solution can be written as 

i & = . 1 ^ F fa-UaA + + q . 1 aA_] (2.103) 

M being the reducible number operator (2.23). Further simplification can be done and 
the expression (2.103) for the new vortex reduces to 

$"7 = f - 1 = , 1 A + + 1 : , 1 A, ) a 2 . (2.104) 

The eigenvalues of the projection operators A± are and 1 and they never contribute 
simultaneously. Owing to this fact (also keep in mind that A± commute with N) the 
expression (2.104) simplifies to 

<$>™ w = - 1 ^^=(A + + A_)a 2 = 1 =a 2 = -^=a 2 (2.105) 

which is same as the n = 2 Witten's vortex. This calculation can be generalized for 
any n and it can be always shown that n-new vortex solution is same as the 2n-Witten 
vortex for all n. It is also easy to show that $oo = , 1 ^=(b^) n is also a solution 
of (2.90) and this solution is same as the fcn-Witten vortex. Note that the leading order 
gauge field and hence the magnetic field are determined uniquely by (A. 18) and (2.88) 
for a given Higgs field. Hence equality of (&™ w , $^f") *° (^oo 5>oo) with vortex number 
kn ensures that the new gauge fields and the magnetic field are also equal to the old 
ones. 
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2.4 Instantons 

2.4.1 Instantons In Commutative Gauge Theories 

Instantons are localized finite action solutions of the classical Euclidean field equations 
of a theory (for a review see [107]). The finite action condition is satisfied only if the 
Lagrangian density of the theory vanishes at boundary. This in turn can lead to different 
topological configurations of the field characterized by its "topological charge" . For Yang- 
Mills theories, the instantons are further classified as Self-Dual (SD) or Anti-Self-Dual 
(ASD) with their topological charges having opposite signs. A simple prescription to 
construct (anti-) self-dual instantons in the Yang-Mills theory is given in [85]. Let us 
first review this construction. We will not distinguish between lower and upper indices 
in this section as the space is Euclidean. 

In order to describe charge k instantons with gauge group U(N) on R 4 one starts 
with the following data: 

1. A pair of complex hermitian vector spaces V = C fe and W = <C N . 

2. The operators B U B 2 e Hom(V, V), I E Hom(W, V), J G Hom(V, W), which must 
obey the equations 

[B u B\} + [B 2 ,Bl\ + // t -J t J = 0, [B u B 2 ]+IJ = 0. (2.106) 



For z = (zi, z 2 ) eC 2 Ri R 4 , define an operator V : V @V @W ->■ V @V as 



2>t 



(yB 2 - z 2 B 1 -z 1 J ) , a 



( -B 1 + z 1 \ 
B 2 - z 2 
J 



(2.107) 
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for anti-self-dual instantons and by 



I , t = ( B 2 - z 2 B x + z 1 I ) , o- 



V 



-B 2 — ^2 
J 



(2.108) 



/ 



for self-dual instantons. Given the matrices B 1 ,B 2 ,I and J obeying all the conditions 
above, the actual instanton solution, A a : W — > W, is determined by the following rather 
explicit formulae: 

A a = ¥d a V, (2.109) 
(a = 1, 2, 3, 4) where ^! : W ^ V @V @W is the normalized mode of the operator 



(2.110) 



Here d a is derivative with respect to the spacetime coordinates x a which are related to 
the ^-coordinates as 



Xi + ix 2 _ X\ — ix 2 x 3 + ix 4 _ x 3 — ix 4 

z l = j z l = 7^ j z 2 — j z 2 - 



V2 ' 



V2 ' 



x/2 ' 



V2 



(2.111) 



For given ADHM data and the zero mode condition (2.110), the following completeness 
relation has to be satisfied 

(2.112) 



X>— i— X> t + ** t = 1. 

It has been shown in [88] that (2.112) can be satisfied even for noncommutative spaces. 
Note that the fields A a are anti-hermitian, consistent with = I. The field strength 
F a p and its dual F a p are given as 



F a p — d a A/3 — dpA a + [A a , Ap], F a p — -e^sF^. 



(2.113) 



The instantons found by the ADHM construction satisfy both the Yang-Mills equation 
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of motion and the (anti-) self-duality condition 



D a F al 3 — 0, F aj3 — ±F a p 



(2.114) 



and has topological charge given by 



Q 



16tt 2 



1 



/ 



d 4 xTr(F a pF al3 ). 



(2.115) 



The trace is over the space W = C . 

2.4.2 Noncommutative Euclidean Space 

To study instantons on a noncommutative R 4 , we will use the notation outlined below. 
The 4-dimensional noncommutative euclidean space is defined by the following noncom- 
mutative coordinates: 



and 9 al3 is a constant anti-symmetric 4x4 matrix. We denote the algebra generated by 
these x a, s by A$. There are three distinct cases one may consider: 

1. 9 has rank (9 a/3 — V a, 0). In this case A# is isomorphic to the algebra of 
functions on the ordinary R 4 . This space may be denoted by R^. 

2. 9 has rank 2. In this case Aq is the algebra of functions on the ordinary R 2 times 
the noncommutative R 2 , which may be denoted by R 2 ^ x R^. Without loss of 
generality, we can choose 



a, = 1,2,3,4, 



(2.116) 







#00 



#000 
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Let us define a system of complex coordinates and a set of operators as 



z\ = -^{x 1 + ix 2 ), z 1 = -Wrr 1 - ix 2 ), z 2 = 4s(x 3 + ix 4 ), z 2 = -4 (x 3 - ix 



V2 
n Zi 



V2 
n t _ z±_ 



V2 V 

a 2 



Z2 



V2^ 

J — n_ 



(2.117) 



which reduces the algebra (2.116) to 

[z 1 ,zi] = 9, [z 2 ,z 2 ] = 0, [a 1 ,a\] = l, [a 2 ,a j 2 ]=0. 



(2.118) 



Here a\ and a\ are like annihilation and creation operators respectively while a 2 , a 2 
are ordinary complex numbers. We can define a number operator by Ni = a\a\. 
The Fock space on which the elements of Ae act, consists of states denoted by 
\rii,z 2 ). Here n\ denotes the eigenvalues of the number operator N x and can take 
only non-negative integral values, while z 2 can be any complex number and denotes 
the eigenvalues of z 2 . 

3. 9 has rank 4. In this case Ae is the noncommutative R 4 . We choose 9 to be of the 
form given by 



g af5 






Q12 













-9 








-9 12 













9 

































-9 








-9 M 













9 






where we have assumed 9 12 = 9 U = —9. Again, we can define a system of complex 
coordinates and a set of operators as in (2.117) but now the algebra (2.116) becomes 



[zi, zi] = [z 2 , z 2 \ = 9, [ai, a\] = [a 2 , a\] = 1. 



(2.119) 



The Fock space on which the elements of Ag act consists of states denoted by 
| ni, n 2 ). Here n\ and n 2 denote the eigenvalues of the number operators Ni = a\a\ 
and N 2 = a\a 2 respectively which can take only non-negative integral values. 
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As already mentioned in section 2.3.1, differentiation in the noncommutative space is 
implimented as an adjoint as in (2. 51), i.e., 

d a f = -i~6^[x^f] (2.120) 

where 9 a p is the inverse of the matrix 6 a p, i.e., 

6 ai 6 1 p = 5 a p. (2.121) 

and is given by 

9 a /3 = —-yflafl (2.122) 

Differentiation with respect to the complex coordinates are 

d z J = ^=[a a J], d- z J = -^=[alf]; a = 1,2. (2.123) 
Again, the integration is implemented by a suitable trace. 

2.4.3 Noncommutative ADHM Construction 

ADHM construction for instantons has been generalized to a noncommutative space 
in [86]. The construction effectively remains same as in the commutative case, only 
change being the replacement of in the right hand side of the first equation of (2.106) 
by the noncommutative parameter 9 for the case of ~R% C x R^ and by 29 for the case of 
R^ c x R^ respectively: 

[B u Bl] + [B 2 ,B%\+rt-jij = 0, [B u B 2 ]+IJ = (2.124) 

for R 2 NC x R^ as in [91] and 

[B u B\] + [B 2 ,Bl] + Ili-fj = 2e, [B u B 2 ]+IJ = (2.125) 
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for R 2 NC x R 2 NC . The Yang-Mills gauge connection A x * is given by 

A xa = -ie a( ,&[x f, ,V] (2.126) 

The gauge field, as in section 2.3.1, is related to the gauge connection by 

D xa = -i6 a/3 xP + A xa (2.127) 



and is given by 

D xa = -i6 a ^x^ (2.128) 

Both A x a and D x <* are again anti-hermitian. In R 2 NC x R^, the components of the gauge 
field along the commutative directions will be given by (2.109), while those along the 
noncommutative axes by D x a. 

Let us first discuss the usual single anti-self-dual U (1) instanton solutions (k — 1, N — 
1) in R 2 NC x R 2 C [90,91] and in R 2 NC x R 2 NC [88-90,95]. For k = N = 1, B 1 , B 2 , 1 and J 
are all complex numbers. As the noncommutative space (described by the coordinates 
z) has translational invariance, we can always choose the origin in such a way that B\ 
and B 2 in r can be taken to be zero. Thus (2.124) or (2.125) ensures that either I or J 
is zero. 1 = gives 

- Z 2 ^ 2 + J£ = 0, -32^1 + 31^ = 0. (2.129) 



while J = gives 

^1-21-02 + ^ = 0, -Z!i/>i- z 2 ^j 2 = 0. (2.130) 



for 



# = 



^2 

V f ) 



(2.131) 
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The two choices are related by a rotation in the plane of complex coordinates z\ — z 2 , 
namely z\ — > —z 2 , z 2 — > Z\. We choose J = without the loss of any generality Now let 
us discuss the two cases seperately. 

R 2 NC x R^: 

Here we get / = \f9 from the equation (2.124). The phase in / does not effect the 
solution for the gauge field and hence has been taken to be zero. The operator (2.107) 
for anti-self-dual instantons becomes 

V = | " 2 ^ ^)=rf( ~ 4 ~' A 1 | (2.132) 
z 1 -z 2 J y ai -a 2 

and its normalized zero mode solution is given by 

1 1 / 8 

ipi = a 2 —=, ip 2 = ai—=, £ = \ — (2.133) 

with 5 = a\ai + a\a 2 and A = 5 + 1. Note that the inverse of the operator A is well- 
defined, but that of 5 is not since \n\ = 0, z 2 = 0) is a zero-mode of 5. 




In this case, (2.125) gives I = \/26 and the operator in (2.107) becomes 

T?t=( v ^ \=V6\ " 4 ~ a{ ^ ]. (2.134) 

ai — a 2 

The zero mode solution is again a 3-element column matrix ^ where we write ip\ = \/6a 2 v 
and i\) 2 = \f0a\v . Then (2.110) becomes 

Av = V29£, v*Av + = ! (2.135) 
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where A = ON = 6{a\ai + a\a 2 ). But this operator does not have an inverse since it has 
a zero-mode 

A|0, 0) = (2.136) 
and hence finding v and £ is a bit tricky. We define a shift operator S such that 

SS* = 1, S^S=1-P, P = |0,0)(0,0|. (2.137) 

Note that although the inverse of A is not defined otherwise, it is well-defined when 
sandwiched between S and . Now we can solve for £ and v: 

£ = $-1^ v = V26-1$-^S\ where $ = 1 + ^ = 1 + 4- (2.138) 

A A N 

which satisfy (2.135). Thus we get 

ii> 1 = V6a 2 V26^r*Si, ij 2 = V9a 1 V29^-^S\ £ = d"^. (2.139) 
We can define the components of the gauge field in terms of the complex coordinates as 



(2.140) 



Di = -j= (f) Xl - iD x ^j , t>i = -j= (d X3 - ib x }j , 

i>i = -j= (b xi + ib X2 ) , b- 2 = ±= (b X3 + w x ^j . (2.i4i) 



Then (2.127) translates to 



D a =\z a + Aa, D- a = - l -z a + A- a) a = 1,2. (2.142) 



with 



A 1 = -^={A Xl -iA X2 ), A 2 = -^=(A X3 -iA x J, (2.143) 
M = (A X1 + iA X2 ) , A 2 = -j= {A x:i + iA Xi ) . (2.144) 
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Also (2.128) translates into 



D = ^ D- 



-Q&Za* = -Dl 



(2.145) 



The solution for k — 1 U(l) ASD instanton becomes 

D a = -l=S$-3a a &Si, D- a = — ]=S&al$-*& 
where the shift operator S, written explicitly, is 



°° / 1 \ 1 

s = i-^K,o>{„;,o|(i- 7 = ai j 

n i — 

fft = l-±(l-a\ |„;,o>(„;,o| 



(2.146) 



(2.147) 



The field strengths are given by 



F af} = -i6 aj [x 7 , A x tt] + i9 01 [a; 7 , A xa ] + [A xa ,A xP ] = -i9 a p + 



D x a , D x f. 



(2.148) 



Their dual are define as 



The ASD condition 



Fa/3 2 ^afiysF^yS 



F a S + F a s = 



(2.149) 



(2.150) 



translates to 



Fn — —F 22 , F\2 — F 12 — 0. 



(2.151) 
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where the field strengths, in our notation of complex coordinates, are given by 



ab 



ab 



Tft [a a , M + -± 



a h-> A a 



F ab - 



± [a a , A b ] - ± [a b , A a ] + [A a , A h ] 
'Te [ a l A b]+7§\ a l A -] +[Aa:A i ] 



ba 



-F. 



ab 



- e Sab + 


D a ,D^, (2.152) 


Da,Db 


(2.153) 


D-a,Db 


(2.154) 




(2.155) 



The equations of motion in the R^c x ^nc m terms of the fields D x a are 

[D xa ,[D xa ,D xP ]] = 0. (2.156) 
The solution (2.146) satisfies both the ASD condition and the equations of motion. 

2.4.4 New Anti-Self-Dual U(l) Instanton 

We can try to get new solutions for noncommutative instantons by using the GBOs. 



Let us first define two operators b and c 



(2.157) 



where Ni is the number operator and A i+ and Ai_ are the projection operators corre- 
sponding to ai- Here b is the generalized operator defined in (2.37). We can easily show 
that 



where N 2 = a} 2 di- 



(2.158) 
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The zero-mode solution \l>o is given by 

V>i = V2c, ij 2 = V2b, £ = 5 ^-±= ai A 1+ + -^== a 1 A 1 ^J . (2.159) 

But this solution is not normalized i.e. ^ Q ^o ^ 1. Usually the single instanton solution 
in R% c x E^, with U(l) gauge group is normalized as [90] 

^ = ^ ^^=- (2.160) 

But in our solution this technique cannot be used because in this case as 

*Stto= ^~_ A j- J(*-l) (2.161) 

vanishes when it operates on the state |0) and the inverse of \T/q^/ does not exist. We 
fix this problem by defining 

V new = y {l-p) i (l-p)u\ wJ = l, v!u=\-p (2.162) 

where p = |0)(0| is a projection operator, and is a shift operator projecting out the 
vacuum. The operator u can be written as 



oo 



u=^2\ni,z 2 )(n 1 + l,Z2\, u 1 — \ni + 1, z 2 )(ni, z 2 \. (2.163) 

ni=0 «i=0 

It should be noted that the new solution in R^ c x is completely non-singular. 



Kc >< ^NC- 



We claim the new solution to be 



D a — ^j^S new Q n £ w b a &newS new i D a — ^—S new ^newb a ^newS new (2.164) 
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where 



$ new = l + — (2.165) 



and 

2 

u m. ■ .u, Yl b i b - ( 2 - 166 ) 

a=l 

Here again the operator is not well-defined otherwise (as M|0,0) = M|0, 1) = 
M|1,0) = M\l,l) = 0), but is well-defined when sandwiched between S new and 
(defined below) which projects out the states |0,0), |1,0), |0, 1) and |1, 1): 

Q Ct - 1 tft C - 1 _ p 

^>new^ n ew ~ x > a new 1 - 'new — 1 1 new ^ 
Pnew = |0,0)(0,0| + |1,0)(1,0| + |0,1)(0,1| + |1,1)(1, 1|. 

The explicit form for the new shift operator is as follows 

= i - EK,o)K,o| (i - -^h) - EW.iXn;,!! (i - 7 ^ T = 1 h) 

71-^ — 71-^ — a 

(2.168) 

The operator b a (corresponding to a a ) is defined as in (2.37). We can check that this 
new solution satisfies the ASD condition (2.151) and the equations of motion (2.156) (for 
details see Appendices A. 5 and A. 6). The topological charge of this new solution can 
be shown to be 4 times the charge of the usual single ASD instanton: Q new = —4 (see 
appendix A. 7). 

This solution is different from the usual k = —4 instanton for U(l) gauge group 
despite the topological charge being the same. We can understand this by observing 
that the difference q between the numbers of a's and a^s in the two solutions is not 
same: q — 1 for the usual ADHM solutions irrespective of its charge and the gauge 
group, whereas q = 2 for the new solution (coming solely because of the operator b in 
the expression of D a given by (2.164)). (The number of a's in \1> is equal to that of a^'s 
in ^ and vice- versa.) The new solution cannot be reduced to the usual instanton by a 
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unitary transformation and hence represents gauge inequivalent configuration. 

If we use the operator b a (z a+ , z a J) defined in (2.41), we can construct an instanton 
solution by repeating the steps we have outlined above. This instanton also has charge 
—4 as the trace of an operator is invariant under unitary transformation. Then the four 
complex parameters z i+ , z 2+ , 2i_, z 2 - can be thought as characterizing the "locations" 
of four instantons with charge —1. It is easy to see that in the coincident limit zi+ = 
z 2 + = Z\- = z 2 - = 0, we recover (2.164). 

We can use the above technique to find a new solution in terms of the GBO b^: 



D = —S ($ )-h( pa) (& ^<?t 

^j-Q*->newY*new ) u a V*new) ^newi 

®ne W = 1 + ^, M = M[ Pl) + M 2 (P2) , M<P«) = b^H^ 
P2-1 oo / 1 \ 

S new = 1 - E KOKil 1 - / ^ » ( 2 - 169 ) 



i=0 n[=0 

with Snew satisfying 



Jm™ + 1 



i= Pl -l,j=p 2 -l 

SnewS new = 1, S new S new = 1 — (2.170) 

i=0J=0 

This solution represents an ASD instanton with charge Q = —p\p 2 . Again, it is gauge 
inequivalent to the k = —pip 2 instanton known in the literature. We could as well have 
used a different shift operator given by 

L = 1"EE \i,n' 2 )(i,n> 2 \ 1 - r -i— 4 P2) • (2-171) 



i=0n' 2 =0 \ a/M 2 (P2) + 1 



Their actions are given by 



; \n±+p u n 2 ) if0<n 2 <p 2 -l 
Siew\ n l, n 2) = \ (2.172) 
\n 1 ,n 2 ) if n 2 > p 2 , 
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S% ew \n u n 2 ) 



n 1 ,n 2 + p 2 ) if < ni < pi — 1 
ni,n 2 ) if ni > p±. 



(2.173) 



Other multi-instanton solutions can be constructed using the reducible representations 
involving squeezed operators (2.47). Our construction of multi-instantonic solutions 
using reducible representations of the standard harmonic oscillator algebra may also be 
generalized for 4k-dimensional instantons as discussed in [108,109]. This exercise will be 
left as a future work. 

In this section we found multi-instantons with charge —p\p 2 (pi,P2 non- negative inte- 
gers) which are not gauge equivalent to known solutions. The charge of the newly found 
multi-instantons has an explicit relation with the representation theory labels pi and p 2 . 
Thus the instanton number does have the information about the reducibility of the space 
along with the topological nature of the solutions. Using the "translated" b operators 
(2.41) we could construct multi-instantons that depend explicitly on pvp 2 number of 
complex parameters. While the full moduli space of noncommutative multi-instantons 
is still not well understood, we hope that this identification contributes partially to this 
question. 



Chapter 3 

Quantum Mechanics with 
time-space noncommutativity 

3.1 Introduction 

In the previous chapter we discussed the theories in noncommutative Moyal spacetime 
where time coordinate remained to be commutative. In such situations where only the 
spatial coordinates do not commute with each other, the quantum theory is conceptu- 
ally straightforward (but nonetheless may display novel phenomena) [110-116]. In this 
chapter we will concentrate on understanding some implications of quantum mechanics 
with time-space noncommutativity, specifically we will work with the Moyal plane E^' 1 . 
We will use the formalism of unitary quantum mechanics on this space as developed by 
Balachandran et. al. [35] (see also [117]). 

When time and space do not commute with each other it is not unreasonable to 
expect that the dynamics of the time dependent processes get altered. We will verify 
this explicitly in the context of a simple model of the forced harmonic oscillator (FHO) 
with the forcing term switched on only for a finite duration of time. In the commutative 
case this is a much studied model. We will compute deviations from the commutative 
case to leading order in 6. These deviations suggest that time-space noncommutativity 
can capture certain nonlinear effects seen in quantum optics. 
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This chapter is organized as follows: In section 3.2 we will briefly review the for- 
mulation of unitary quantum mechanics on R^' 1 [35]. In section 3.3 we will solve the 
problem of the FHO perturbatively in 9 and compute corrections to the transition prob- 
abilities between simple harmonic oscillator (SHO) states. These corrections suggest the 
noncoherent nature of the time-evolved vacuum state and are the reminiscent of those 
seen in nonlinear quantum optics [118]. To flesh out this analogy better we study the 
time-evolution of uncertainties in position and momentum in section 3.4. Encouraged 
by these results we, in section 3.5, suggest a correspondence between the nonlinearity in 
quantum optics and the quantum mechanics on R^' 1 . 

3.2 Unitary Quantum Mechanics on R^' 1 

The noncommutative space R^' 1 is described by the coordinates x^s satisfying 

[x^,x u ] = iQe^y with e^ u = -e VfM and e 01 = 1, (3.1) 

where \i and v can take values 0,1. Without loss of generality we can take 9 > 0, 
as its sign can always be flipped by changing x\ to —x±. Let ^(R 1,1 ) be the unital 
algebra generated by x and x±. We associate to each a G Ae (R 1 ' 1 ), its left and right 
representations a L and a R : 

& L f3 = ap , & R f3 = f3a , $ e A e (R 1,1 ) . (3.2) 

Unless stated, we work with the left representation. 

For a quantum theory, what we need are: (1) a suitable inner product on Ae (R 1 ' 1 ); 
(2) a Schrodinger constraint on Ae (R 1,1 ); and (3) a self-adjoint (with respect to the inner 
product defined) Hamiltonian H and observables which act on the constrained subspace 
of A e (R 1 ' 1 ). 

1. The Inner Product: 
There are several suitable inner products and they are all equivalent to each other as 
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discussed in [35]. Here we discuss only one such example. Using the commutator (3.1) 
any a G ^(R 1,1 ) can always be written as a = J d 2 k a(k)e tklXl e lkoXo . We associate a 
symbol 0:5 corresponding to each such a given by 

a s (x ,x 1 ) = J d 2 k&(k)e tklXl e ik ° x ° . (3.3) 

Note that x and x\ and hence 0:5 are purely commutative. The inner product is defined 

as 

(a, /#) t = J dx 1 a* s (t,x 1 )(3 s (t,x 1 ) ■ (3.4) 

2. The Schrddinger Constraint and time evolution: 
The operators p and pi, given by 



— = A --ad -i— 
8xq ^° 9 X1J dxi i " 1 9 



i ^Z~ =Po = -Tjadxu -*-HT" = Pi = -^ ad £o) ( 3 - 5 ) 



generate time and space translations respectively. The adjoint action is defined as 

ad^ = [i,# (3.6) 

It can be easily checked that the canonical commutation relations are satisfied: 

[x^Pv] = -iri^ with 77^ = r] Ufl and 7701 = 0, r] 00 = l,r] n = -1 (3.7) 

The Hamiltonian H, in general, may depend on x^, Xq and p\. The possible dependence 
of £f and Xq can be bypassed by 

Xi = 9p + x\ , x% = 8p! + x§. (3.8) 

Also, there is no dependence on p Q assumed in the line of the commutative case where 
there is never such dependence of H on id Xo . Now note that the inner product (3.4) 
has an explicit dependence on the parameter t and hence there exist more than one null 
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vectors with respect to this inner product (actually any vector which vanishes at xo — t 
is a null vector). But this fact need not bother us as we are only interested in those 
states that satisfy the Schrodinger constraint 

(po - #) i> = 0. (3.9) 

The Hamiltonian H depends on Xi,Xq and p\. Since Xq commutes with x\ and pi we 
will choose x§ as "time" . 

It is easy to write down the formal solution of the Schrodinger constraint and find 
the time evolution. The time evolution is given by xq — > Xq + r (or equivalently by 
Xq — > Xq + r) . Thus the amount of time-translation is always commutative, though the 
time-operator itself is noncommutative. The time evolved wave functions satisfying the 
Schrodinger constraint are of the form ip(x ,Xi) = U (£^,77) x(xi), where 

(3.10) 

As Xq commutes with x\ and pi, the time dependence of the time evolved wave function 
given above will mimic the same for the commutative case 9 = 0. We know in the 
commutative case a wave function which vanishes at some finite time t, will vanish for 
all times. Hence the only null vector satisfying the Schrodinger constraint for the case 
of nonzero 9, is tp = and there are no other non-trivial null vectors! 
3. The Spectral Map: 

Consider a time- independent Hamiltonian H — + V(xi). The corresponding commu- 
tative Hamiltonian is H = — ^^2+V( x i), with eigenfunctions ip E (x ,Xi) = ifE{xi)e~ lE:> 
and eigenvalues E. The spectrum of the corresponding noncommutative H will be given 
by tpE = e~ tEx oif E [xi) = (pE(xi)e~ tEx ° with the same eigenvalues E as H<pe(xi) = 
Eip E (xi). Here (Pe(&i) has been obtained by replacing x 1 with x\ in ip E ( x i)- 



U(x§,t t ) = Texp 



— 1 



/ drH(r,x L , pi) 

Jtj 
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3.3 QFHO in and their Transition Probabilities 

Let us recall the dynamics of a QFHO in ordinary spacetime. For a detailed discussion 
one can look into the section 14.6 of [119]. The Hamiltonian of this system is given by 

H(t) = ^- + l -m Q u 2 x 2 + f(t)x + g(t)p, (3.11) 

where mo is the mass of the particle and u is the angular frequency of the oscillator. We 
are interested in real functions obeying 

f(t),g(t) = 0foit->±oo. (3.12) 

At t — > — oo the Hamiltonian is simple harmonic and we assume the system to be in one 
of the eigenstates of this SHO Hamiltonian. At t — > oo the Hamiltonian again becomes 
simple harmonic and we try to find the probability (the transition probability) for the 
system to be in any arbitrary eigenstate of the SHO Hamiltonian subjected to the fact 
that the system was in some already given eigenstate at t — > —oo. For this what we do 
is the following: 

• First we assume our system to be in an eigenstate (f> n (x) at t — ti — > — oo. 

• The state (f> n (x) evolves under the SHO Hamiltonian from t — ti — > — oo to t — Ti. 

• At t — Ti the interaction gets switched on. 

• The system then evolves under the full Hamiltonian (3.11) from t — Ti to t — T 2 . 

• At t = T 2 the interaction gets switched off. 

• The system again evolves under the SHO Hamiltonian from t — T 2 to t — tf — > oo. 

• We find the inner product of the final state we get at t = tf — > oo with the 
eigenstate (f> m (x). This gives the Transition Amplitude A mn while its absolute 
square gives the Transition Probability P mn - 
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The generalization of the above Hamiltonian in R^' 1 is 

p 2 1 1 - 

H = + -m Q u 2 x\ + - [f(x )x 1 + x 1 f(x )] + g{x )pi = H + E u (3.13) 



with 



p 2 1 1 
Hq = + -m w 2 ^, i7j = -[/(x )xi + x 1 f(x )] + 5-(£ )Pi- (3.14) 



As x and commute with each other, the ordering does not matter in the last term. 

To define the transitions for the above Hamiltonian consider the time evolution by 
an amount r. The functions f(xo) and g(xo) have the properties of vanishing in the far 
past and the far future, i.e., 

/, g(x + r) — > as r — > ±oo. (3.15) 

We shall find the transition probabilities (P m ,n) for an SHO state "n" at initial time 
(r — > — oo) to go to some other SHO state "m" at final time (r — >■ +oo) after evolving 
under the Hamiltonian (3.13). The Spectral Map tells us that the energy spectrum of 
the SHO Hamiltonian in R^' 1 is same as that of the commutative one, i.e., 

E n = hw(n+^j, M^x l ) = <p n (x 1 )e-^ n+1 *>\ ( 3 . i6 ) 

where 4> n {xi) is the eigenfunctions of the commutative SHO Hamiltonian. The orthon- 
armality of the eigenfunctions ^ n {xQ,X\) with respect to the inner product defined in 
section 3.2 can easily be checked and is shown explicitly in Appendix B.l. 

The transition probabilities for our problem can be found by computing the same for 
the commutative Hamiltonian obtained after replacing 

x x ^x, pi->p, x = Xq = --pi + x$ -> ~p + t, (3.17) 

in the Hamiltonian (3.13). Here t has come in place of the "time" Xq which commutes 
with xi and p±. Also, we have retained the fundamental constant h while it was taken to 
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be 1 in section 3.2. To linear order in 6 we obtain the following commutative Hamiltonian 

H(t) = H + Hj(t) = H + H I0 (t) + 6H n (t), (3.18) 



with 



H — mr)U 2 x 2 = hou { a} a H — ) , (3.19) 
2m 2 ° V 2 J 

H I0 = f(t)x + g(t)p = z*(t)a + z(t)a\ (3.20) 



and 



Hn = ^(-g'(t)p 2 - ^f'(t)(xp + px) 

i 

h 



-"' 1 z (t)a — z \t)a) + i\ — - — g (t)(2a'a + 1) (3.21; 



The function z(t) is related to fit) and g(t) as 



*(*) = ^ (/(*) + imaug(t)) . (3.22) 
Also, a and a* are the annihilation and creation operators respectively defined as 

=S> (3.23) 

The nonlinearity in the Hamiltonian (3.18) is purely due to the noncommutativity. This 
provokes us to model certain types of nonlinear phenomena in quantum optics by the 
noncommutativity between time and space coordinates. This analogy will be further 
studied in section 3.5. Let us now continue with calculating the transition amplitude 
which is given by 

Am, n (tf> T X T uU) = {KWlitf^Uoitf^U^T^U^U)^), (3.24) 
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where U (f, t) and U (t' , t) are the time evolution operators from time t to time t' for the 
Hamiltonians H and H{t) respectively, i.e., 

U (t',t) = e -* Ho ^-*), U{t',t) = T [ e -s/t*'*/f(r)j ? (3>25) 

the latter one being the time-ordered exponential. This gives 

An,^*/,^;^^) = et^^-^+^^-^)](0 m |f/(T 2 ,T 1 )|0 n ). (3.26) 

The state \tp(t)) = U(t,Ti)\<f) n ) evolves according to the Schrodinger equation for the 
Hamiltonian (3.18) 

^j t ~H^j \m) = H I (tM(t)), (3-27) 

with the initial condition \ip(t = Ti)) = |0„). If we define the Green's operator function 
G(t,t ) as 

ih%- - H ) G(t, t ) = 5{t - to), (3.28) 



dt 

then solution of the Schrodinger equation (3.27) will be 



\m) = 10(0) + / dt G(t,t )Hj(t )m )), (3.29) 
which in turn gives the Born series 



+oo 



m)) = 10(0)+/ dtoG^toWtoMto)) 

J — oo 

/+oo p+oc 

dt dtiG(t,to) J ffj(t )G(to,ti)^/(ti)|^i)) + .- (3.30) 
■oo J —oo 

Here | </>(£)) is the solution of the homogeneous equation 

ihj t - H \ |0(O) = (3.31) 
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which is nothing but the Schrodinger equation for SHO. G has been found in the Ap- 
pendix B.2 (see (B.8)). Note that the 6-function in the expression of the G restricts the 
integration over tj in (3.30) within the limit of — oo to tj-i (i_i = t). Thus, at t — Ti the 
integrations are only in the intervals when the interaction was switched off, i.e., Hi = 0. 
Hence, we get 

|^(t = T 1 )) = |0(t = T 1 )) = |0„). (3.32) 

The solution of the homogeneous equation (3.31) with the initial condition (3.32) is given 
by 

|0(t))= e -^(^)|0 n ). (3.33) 
Putting this back in (3.30) we get 

/+oo 
dt G(t,t )Hi(to)\(f>(to)) 
■oo 

/+oo p+oc 
dt dhGit^H^G^tjHjihMh)) + .... (3.34) 
-oo J —oo 

For t = T 2 we get 

/+oo 
dtoGiT^toWtoMto)) 
-oo 

/+oo r+oo 
dt / dtiG(T 2 ,i )^/(t )G(*o,*i)ffj(ti)|^i)) 
-oo J —oo 

+.... 

(3.35) 



We put this back in (3.26) to get 



A m , n (tf, T 2 ; Ti, ti) = ^2 B j( t f' T 2; 7i, U), (3.36) 

j=0 



with 



Bo(tf,T 2 ;Ti,ti) — 5 rnn 



(3.37) 
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and 



/+oo p+co r+oo 

dt / dti... / dt j -iF^ n (tf,T 2 ;t ,ti,...,tj- 1 ;Ti,t i ) 
-oo J— oo J — oo 

(3.38) 

for j = 1, 2, .... Here 

Fi hn (t f ,T 2 ;t ,t 1 ,...,t j . 1 - t T 1 ,U) = eK^^-*') + ^^- Tl >](^ B |G(r 2 ,to)fr / (to) 

G(t o ,ti)^(ii)...G'(^_ 2 ,^_ 1 )// 7 (^_ 1 )|0(^_ 1 )) 

(3.39) 

We put the expressions for all G(t — i')'s as found in (B.8) to get 

Fi >n (t f ,T 2 ;t ,t 1 ,...,t j - 1 ;T 1 ,t i ) = (-^ 6(T 2 - i )e(f - *i)...efo_ 2 - f,_i) 

^(^Oe*^- 1 !^--!)). (3.40) 

The if mt 's are defined as 

tf™*(t) = e^H^.^e-^. (3.41) 
We put 10(^-1)) = e-s^*fe-i- T O|0 n ) ( see (3 33)) to get 

FL,n{tf-> T2] to, ti, t-,_i; Ti, tj) = e(T 2 -*o)e(*o-ti)...e(* J -_ 2 - 

e^^-^^^l^r'^o)^^!)-^^-!)!^)- 

(3.42) 

Writing Hi{t) = H I0 (t) + 6H n (t) (see (3.18)) and hence H} nt (t) = Hj^(t) + 0iZ}?*(t) to 
separate out the ^-dependent and independent parts up to linear order in 0, gives 

A m ^t h T 2 -T x X) = e^-^i^A^iT^T,) + 0AW(T 2 ,T 1 )}\<f> n ) , (3.43) 
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with 



/>H) 



1 / dt [--)Q(T 2 -t )H^(t ) 



-oo 

/+oo /.+oo / ^ \ : 



-oo ■> — x 

+ ... 



/ d*i ( J 6(T 2 - t )e(t - h)H% \t )HX (h) 



= T 



e ft J -°° 



(3.44) 



• /-+00 

^(T^TO] = --J dt G(T 2 -t )Hy*(t ) 

-z) y *o y dti6(T 2 - * )e(*o - *i) 



.[if}? (t )^lS*(*i) + H^it^H^it,) 

3 /-H-oo /"+00 /"+oo 



(2 \ /"+°° /-+oo /-+oo 

--) / dt dh dt 2 e(T 2 - to)e(t - *i)e(ti - t 2 ) 
'"J J —oo J —oo J —oo 

. [H}f (foWfaWte) + fliS*(*o)fll? (*i)fl}?(* 2 ) 

+^(to)^(ti)^f(t 2 )] 
+.... (3.45) 

Hi (t) and Hn(t) are again defined in accordance to (3.41). The arguments of A(°)(T 2 , Ti) 
do not come in the formal expression as the integrands of all the integrals vanish for 
t > T 2 . The expression for A^(T 2 ,Ti) can be simplified to (see Appendix B.3) 

i f 00 _i 

A (1) (T 2 ,Ti) = -lAW^.Tx) y rft [AW^TO] " 1 ^(io)A(°)(t ,T 1 ). (3.46) 
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where A^(t,t') with arbitrary arguments is defined in (B.ll). Putting this in equation 
(3.43) we get 



Am,n(tf,T 2 ;T 1 ,t i ) = en 



(E n — Em)ti 



l-j:0 I (fto[A(°)(t ,r 1 )]- 1 fllf(to)A( )(to,r 1 ) 



\4>n). 



(3.47) 



Here we have removed the arguments of A^(T2,Ti) as they don't come in it's formal 
expression. For the next few steps we are going to use the following identity extensively: 



e xA Be -xA = B+±[A,B] + ^ [A, [A, B]] + ^ [A, [A, [A, B]}] + ... 
First of all to find H l ^{t) and Hf^(t) we calculate the following commutators: 



(3.48) 



[H ,H I0 (t)] 
[Ho, [Hq, H I0 ( 



[H ,H n (t)] 
[H ,[H ,H n (t)}} 
[Ho,[Ho,[Ho,Hn(t)]]} 



hu(-z*(t)a + z(t)a!) (3.49) 

h 2 u 2 H I0 (t) (3.50) 
...and so on. 

-2tu^^ {z*'(t)a 2 + z'(t)a,V) (3.51) 

ifiw 2 (H n (t)+m g'(t)H ) (3.52) 

4hw 2 [Ho,H n (t)] (3.53) 
...and so on. 



Now, using above commutators in the identity (3.48) gives 



H%\t) = e* Hot H IO (t)e-* Hot 

= z*(t)e- iuit a + z(t)e iu ' t a ! 



(3.54) 
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and 



Hjf(t) = e^ Hot H n {t)e-i Hat 



m Q hw 

h\l 2 

m co , 



z*'(t)e 



-2iu)t 2 



g'(t) (2a f a + l) 



(3.55) 



We put the expression of Hf^t) given above to find A^(t : to) (see (B.ll)) and follow 
the discussions given in pages 338-340 of [119] to eliminate the time ordering. This gives 



^(°)(^ ) = e ^(t,to) e -e(t,to)aH(t,to)a^ 



(3.56) 



with 



0(Mo) = ^ f dr^ £dr 2 [z* (r^z^e^^ - z^z* (r 2 )e^-^} (3.57) 



and 



f(Mo) = -fc / dTz(r)e 



h 



(3.58) 



to 



The limit t - >■ —oo,t — >■ oo gives the expression for A^ ). Now to simplify the integrand 
in (3.47) we find the following commutators 

[r^cTOa-^cTOat,//}-*^)] 



2 / mpfru) 

h\J 2 



2 



g'ito^ito,^) +iz*'(t )Z(t ,T 1 )e—» | a 



-2iajto \ 
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and 

[C(t , T,)a - £(* , TOat, [£*(i , T,)a - £(t , T^, H%*(t )]] 
-2rrW(*o)|£(*o,:Ti)| 2 + iV2^^ 



All higher order commutators vanish! Thus we get 

+a 2 (t ,T 1 )a + a;(t ,T 1 )a t 

+a 3 (t , Ti)a 2 + 0:3 (t , Ti)a t2 

+0:4(^0, Ti)a^a (3.59) 



with 



ai(*o,7i) = ^-^(*o)-mow^'(to)|e(*o,7i) 



+ y ^ ^'(to)e 2 (to,T 1 )e- 2 ^ - ^ / (t )r 2 (to,T 1 )e 2 ^] (3.60) 
"2(^0,^1) = -m w^ / (*o)r(*o,Ti) 

+j)/^z"(to)Z(to,T 1 )e-*^ (3.61) 



asfo.TO = y^z'Me-™* (3.62) 
"4(^0,^1) = -m ug'(t ) (3.63) 

Finally we get the expression of the transition amplitude as 

+/3 2 * V / ^Tl J D m ,n+l(0 + /W™(™ - l)An,n- 2 (0 

+/3 3 V(n+l)(n + 2)D m , n+2 (0}] , (3.64) 
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with 



= ~T dre lWT z(T) (3.65) 

J-oo 
■ roc roc 

= 2VJ dTl J rfr 2^(n)^(r 2 )e-^ (ri - r2) -^(r 1 )^(r 2 )e^-^)] 

= real (3.66) 



>\ — 



-m u l^drg'(r)\ar)\ 2 + 1 -^ 



i m huj 



2 

-+oo 



2 



/+oo 
dr [z*'[r)i 2 {T)e- 2wr - z '(r)f 5 \r)e 2wr \ (3.67) 
■00 

-m G u £J drg'(r)C(r) + f /^/^J dr z^r)^' 2 ^ (3.68) 



-00 

-+oo 



A = I drz*>(r)e— (3 . 69 ) 



ft V 2 

Here the function is given as 

Z(t) = -±f t dre^z(r). (3.70) 

An,n(£)' s are the matrix elements of the displacement operator £)(£) = e -£*<H-£a t gj ven 
by [120] 

A»,»(o = \/4 e ^ l5|2 r- n ^r n (iei 2 ), (3.7i) 



ni _i| C | 2 

— f 
m! 



L„(x) are the associated Laguerre polynomials. Also, the limits of the integrations have 
been extended to — oo and oo as the integrands are zero in the extended region. Note 
that the contribution from a 4 given in (3.63) to the transition amplitude vanishes by 
virtue of (3.12). The transition probability is given by 

Pm,n = \^m,n(tf,T 2 ;Ti,ti)\ 2 (3.72) 



as usual. The arguments have been omitted as the formal expression of transition prob- 
ability does not contain the times t/,T 2 ;Ti,tj. 
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3.3.1 n = 



For vacuum |0 O ) as the initial state, the transition amplitude is 



A mfi (t f ,T 2] T u U) = e^e^-^e * 



'ml 



i-jem-PlC + PlC 2 ) 



+ m J^2^C - 2/3 3 T 2 ) + m(m - l)^/^* 2 } 



(3.73) 



We have used the explicit forms of the following associated Laguerre polynomials: 



L™(x) 
x) 



= 1 



Lm—2 
2 



m — x 
1 



\x 2 — 2mx + m[m — 1)] 



(3.74) 
(3.75) 
(3.76) 



The transition probability becomes (upto linear order in 6) 



P-mfi — \A m fl(tf,T2]Ti,ti 



i2 = e -\e\^l 

ml 



1 + + mA 2 + m(m - 1)A 3 } 



, (3.77) 



with 



A^ImifoO-ImiM 2 ), A 2 = ^(2/ m (^ 2 )-M«)), A 3 = -^Im(M 2 )- 

(3.78) 

Note that as m — > oo, the ^-correction starts dominating and in this case the expansion 
upto linear order in 9 is no more meaningful! Hence, the above result is valid only for 
those m- values which are far smaller than l/y/m u8 (in the unit h = 1). For 9 — > 0, the 
transition probability becomes the well known Poisson distribution as expected. 

As a specific example let us work with the functions fit) and g(t) of the form (see 
Figure 3.1) 

/(f) = f [&(t + T)-&(t-T)] 



git) = go [e(t + T)-e(t-T)} 



;T > 0. 



(3.79) 
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This gives 



z (t) = zo[e(t + T)-e(t-T)], z 



2m u 



[fo + imawgo] (3.80) 



Their derivatives are 



x'(t) = x [6(t + T) - S(t - T)] , iorx = f, g, z. 



(3.81) 



For these functions we calculate to be 



m = 



_ 2izp sinuT 
fru) 



t < -T 
—T <t<T 
t > T 



(3.82) 



We again calculate the following quantities: 



?2 = t: sin coT 
a 



im g z 



U \2-2iuiT 

~T — \ z 0\ e 

nw 



(3.83) 



2moW * . 
P3 — —\l ^ z sm2wT 

2izo sin ojT 



which further give 



(3.84) 
(3.85) 



| 1 2 sin 2 uT [Re(zo) cos 2uT + Im(z ) sin 2uT] (3.86) 



ffioj V hcu 



and 



Im(M 2 ) = 



4 2m 



° l z \ 2 sin 2 uoT sin 2uT Im(z ) 



h 2 co V ^cu 



(3.87) 
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Thus we get 



A 1 = 2 {° 5 (/o + ^n^ 2 ^n) sin 2 cos 2uT (3.88) 
m nuj 6 



A 2 = m Logo sin 2wT — f cos 2o;T (3.89) 
/o + ^o w2 ^o 



A 3 = - - 2 ,° cotcT. (3.90) 



The choices m = 1, u; = 1, /o = V^, go — v^5 and T = | (in natural units, i.e., ft = 1) in 
commutative case (0 = 0) give the following Poisson distribution: P mfi = e~ 20 ^, while 
for nonzero 9 the probability distribution modifies to P m fl = e~ 20 ^p[l + 2^v / 5(m — 20)]. 
The ^-correction becomes of the order of the ^-independent part when m approaches the 
value m(9) = ^20 + ■ Hence, our result is valid only in the region where m < m{9). 
Note that m{9) ~ | rather than ^ because the A 3 is identically zero for the choices 
taken. We choose 9 = 0.01 (m(9 = 0.01) « 42) and get 



m,0 



20 "" [1 + 0.02^^-20)]. (3.91) 



mi 



This deformed distribution along with the Poisson distribution is shown in Figure 3.2. 
Such deformation of the Poisson distribution suggests that the vacuum does not evolve 
to be a coherent state anymore. To explore this further let us look at the time-evolution 
of position and momentum uncertainties. 

3.4 The time evolution of A x and A p 

The expectation value of any operator O in a state ip(x , x±) at any time t is defined to 
be 

(6} t = ($M) t - (3-92) 

Also, 

(0) t+T = (j>,6ty = ^(x + r,xi),(9^(xo + r,xi)) t . (3.93) 
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Figure 3.1: The behaviour of functions f(t) and g(t) with t 

Figure 3.2: The modified distribution (8 = 0.01) along with the Poisson distribution 
(9 = 0) for m = 1, u — 1, fo — \/5, go = \/5, T — | and 7i = 1 

Hence the time evolution of the expectation value of an operator is given by that of the 
state in which it is being calculated. For the QFHO in Rg' 1 the time evolution of any 
operator O will be given by 

d (6) = ^(d) + U\H(t),d\), (3.94) 



dt dt h 

where H(t) is the Hamiltonian (3.18). The uncertainty in any observable O is given by 



A a = J(0*)-(Oy. (3.95) 



For an observable not having an explicit dependence on time we will have 

jtf = % - (([H(t),d 2 ]} -2(d)([H(t),d])) (3-96) 
To find the evolution of and we calculate the following commutators 

[H(t),x] = % [29g'(t) - — J p + i6f{t)x - ihg(t) (3.97) 
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[H(t), x 2 ] = i \26g\t) - — j (xp + px) + 2i9f'(t)x 2 - 2ihg(t)x 
[H(t),p] = ihm u 2 x - i9f(t)p + ihf(t) 
[H(t),p 2 ] = ihm u 2 (xp + px) - 2i6f'{t)p 2 + 2ihf(t)p 



(3.98) 
(3.99) 
(3.100) 



Thus the evolution of A 2 , and A 2 is 



dt^x 



1 20 // 



n g>(t))(l(xp + px)-(x)(p))-ff>(t)A 
A 2 p = -2m u; 2 (l(xp + px)-(x){p)) + ff(t)A 2 p 



iA2 

dt 



(3.101) 



Let us define 



Using 



A xp = -(xp + px) - (x)(p), 



(3.102) 



[H(t),xp + px] = 2% \2Qg\t) - — j p 2 + 2ihm u 2 x 2 + 2ihf(t)x - 2ihg(t)p (3.103) 



we find the following first order coupled equations 



^A 2 

dt^x 

dt^P 
A A 

dt^ x v 



2 - f ✓(*)) - f /'(*)* 

-2moo; 2 A^ + f f'(t)A 2 p 
k- 2 i9'(t))A 2 p -m u 2 A 2 x 



(3.104) 



For vacuum as the initial state, the initial conditions for the above are 

A 2 x (t -oo) = 7 ^—, A 2 v (t -oo) = 7 ^-, A xp (t -+ -oo) = 0. cim.; 



2m w 



Our strategy for solving these equations is simple. We do so perturbatively in 6. We 
write the perturbative expansion of the solution up to linear order in 6 



Al = Af) + 0Af >, A 2 p = Af ) + 9Af\ A xp = Ag + 0A« (3.106) 
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Splitting the equations in the part independent of 6 and the part which is linear in 6 we 

get 

Ul {0) = (3-107) 



dt 

-A<°) = 

with initial conditions 



dt x m 

^Af) = -2m c 2 A(° p ) (3.108) 
jUg = -Af ) - m , 2 Af ) (3.109) 



Af)(t^-oo) = ^_, Af)(t^-oo) = ^, A£>(f->-oo) = (3.110) 



and 



^-Af) = — A« - ^(t)A£) - |/'(*)Af ) (3.111) 
|Af) = -WAg + |/'(t)Af > (3.112) 
|A« = i-Af)-^'(t)Af)-m o; 2 Af) (3.113) 

with initial conditions 

Af)(t^-oo) = 0, Af)(t^-oo) = 0, Ag>(f->-oo) = 0. (3.114) 

The solutions of the 6 independent part with the given initial conditions can be easily 
found to be 

A ? 0) (0 = 2^J, Af)(t) = ^, Ag(f) = (3.115) 

for all time t. The results are also obvious from the fact that for the commutative case 
(9 = 0) the time-evolved vacuum state is nothing but a coherent state in which the 
uncertainties are constants and are given by the above values. Putting above in the 
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equations for the #-part, we get 



_l A 2(i) 
dt x 

_1 A 2(1) 
dt P 

dt xp 



The first two equations give 



= -AS-— fit) 
-- -2m w 2 Ag + m uf'(t) 



m 



-A 2 « - m ujg'(t) - m cu 2 A 2 x V 



dt 



whose solution satisfying the initial conditions is given by 



(3.116) 
(3.117) 
(3.118) 



(3.119) 



Af ) = -m^Af > 



(3.120) 



We put this in (3.118) to get 



d 
dt 



Ag = -2m u 2 A^ > - m Q ug\t) 



(3.121) 



Differentiating (3.116) once and using the above equation we get the decoupled equation 
in A 2 / 1 ^ as 

(li* + A ' W = "WW - —f"(t) (3.122) 
\dt z J m uj 

with the initial conditions 



Af)(^-oc) = 0, ^Af)(t) 



= 



(3.123) 



t—t— oo 



Here, the second condition comes from the equation (3.116) and the initial conditions of 
Axp(i) and /'(£). The solution is 

A 2 «(t) = —f(t)-—[ drsm{2uj(T -t)}f(r) - 2u [ dr cos{2w(r - t)}g(r) 

(3.124) 
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Using above in (3.120) and (3.116) respectively we get 



AjW(i) = m uj f(t) + 2uf drsm{2uj(T -t)}f(r) 

J — oo 



9(r) 



m LU 



+2m u 2 / dr cos{2w(r — t)} 

J —oo 

I drcos{2uj(T -t)}f(r) 

J — oo 

2u I dr sin{2w(r - t)}g(r) - g(t) 

J —oo 



1 '" 

moJ-oo 



Thus, the uncertainties upto linear order in 9 are 



(3.125) 



(3.126) 



h 9 itlq 



2m u 2\2hu 



2 

m 



a rt 

f(t) + — dT S m{2io(r-t)}f(r) 

m J-oo 



+4u 2 j dr cos{2w(r - t)}g(r) 

J —oo 



(3.127) 



m hu 9m /m uj 



2 ' 2 V 2h 

-t 



2 

TO 



f(t) + — f t drsm{2u(r-t)}f(r) 

"W-oo 



+4u 2 ( dr cos{2u(t -t)}g(r) 

J —oo 



(3.128) 



A xp (t) = 9m u 



-g(t) + — / dr cos{2w(r - t)}f(r) - 2u [ dr sin{2w(r - t)}g(r 

m 0j-oo J-oo 

(3.129) 

The fundamental uncertainty product (to linear order in 9) is 



A x {t).A p {t) 



h 



(3.130) 



Thus the vacuum state evolves to a "squeezed state" rather than a coherent state as in 
the commutative case [121]. The uncertainties in the commutative case depend only on 
the product m u. But, their ^-corrections change with oj even if mou is kept constant. 
Also, the squeezing effect is oscillatory in time as is obvious from the ^-dependent terms 
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in (3.127), (3.128) and (3.129). For the specific forms of fit) and g(t) of (3.79) we get 



A x (t) = { 



; t<-T 



+2m Q g Q u sin{2w(t + T)}] 

+ 2m g u (sin{2w(t + T)} - sin{2w(t - T)})] 



-T <t <T 



; t>T 



(3.131) 



A p (t) = < 



and 



2 

rnphu) | /mow 



+ fv / W[ 2 /ocos{2a;(t + T)} 
+2mo5'oW sin{2w(t + T)}] 

+ fv^[2/o (cos{2a;(t + T)} - cos{2o;(t - T)}) 
+ 2m g uj (sin{2w(t + T)} - sin{2w(t - T)})] 



t < — T 



-T <t <T 



t > T 



(3.132) 



A xp (t) = < 



■ t<-T 

|[2/ sin{2a;(t + T)} 

-2m ug cos{2w(t + T)}] ; -T < t < T 

| [2/ (sin{2w(t + T)} - sin{2w(t - T)}) 

-2m ujg (cos{2w(t + T)} - cos{2w(t - T)})] ; t > T. 



(3.133) 

For the choice of parameters m = 1, u = 1, fo = \/b, go = ^V^, T — | and = 0.01 in 
natural units (ft = 1) we get 



A x (t) 



i 

^ + 0.01 
^ + 0.01 



• t < — - 

5 (cos2t + |sin2i) ;-§<£<§ 
|sin2t ;£>!, 



(3.134) 
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A p (t) 



i 

V2 



t < 



j= - 0.01 */§ (cos2t + | sin2t) 
^- 0.01^1 sin 2t ;f>f 



and A xp (t) = < 



O.Olv^ Qcos 2t-sin 2t) 
-0.02 v / 5sin2t ;t>f. 



< t < I 



;t< 



-!<*<! 



(3.135) 



(3.136) 



Figures 3.3 and 3.4 show the time-dependence of the different uncertainties. The discon- 
tinuities at t — ±| is simply the manifestation of the fact that the functions f(t) and 
g(t) themselves are discontinuous at these times. Before the interaction was switched 
on, the uncertainties were having values equal to those for the vacuum state. During the 
time of nonvanishing interaction (and even after the interaction gets switched off!), they 
oscillate with frequency equal to twice that of the oscillator. 




2n ?» l n 




£ 2n 5£ in 



Figure 3.3: The time-dependences of the uncertainties A x and A p for m = 1, oj 
/„ = V5, g = |V5, T= f, 9 = 0.01 and H= 1 



Figure 3.4: The time-dependence of A xp for the same choice of values 
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3.5 Implications in Quantum Optics 

In Quantum Optics a monochromatic (single-mode) coherent light field is usually de- 
scribed by the harmonic oscillator coherent states [122]. It has also been shown that a 
coherent state (in particular the vacuum state) remains to be coherent under the FHO 
Hamiltonian [123]. The annihilation and creation operators for photons are related to 
the field quadratures X\ and X 2 by 

a = X x + iX 2 , a f = X 1 - iX 2 , (3.137) 

Xi and X 2 being hermitian. The commutation [a, a 1 "] = 1 translates to [Xi,X 2 ] = \. 
The coherent state has different uncertainties as A Xl = |, A X2 = \ and A Xl x 2 — =^ 
Ax 1 -Ax 2 — \ which is the minimum. Also, the photon count (probability for having a 
certain number of photons) in the coherent state is given by the transition probabilities 
of the corresponding number eigenstate and the profile is Poissonian. 
The FHO Hamiltonian 

/ oh 

H(t) = hu(X* + X%) + J f(t)X! + v / 2hm tug(t)X 2 

V m uj 

= ftw ^aW +z*(t)a + z(t)a ( (3.138) 

(z(t) is related with f(t) and g(t) by (3.22)) with the effective noncommutativity between 
time and the field quadrature X± of the form 

[t,X 1 ]=i ) J^e (3.139) 

will allow us to use the calculation of the previous sections. The photon count will be 
given by (3.77), while the uncertainties in the field quadratures = \j {Xfj — (Xi) 2 , 
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A X2 = y/(X%) - (X 2 ) 2 and A XlX2 = + X 2 X 1 ) - (X 1 )(X 2 ) will get modified as 



A*(f) 



1 Qmr 



2 Ah 



f(t) + — f drsm{2io(T-t)}f(T) 

m 0j-oo 

+4oo 2 f dr cos{2w(r - t)}g(r) 

J — oo 



2 

m 



(3.140) 



Ax 2 (t) 



1 9m 
- H 

2 Ah 



2 

m 



f(t) + — f t dr S m{2u(r-t)}f(r) 
+4w 2 / rfr cos{2w(r - t)}g(r) 

J — oo 



(3.141) 



AW*) 



9m u! 
2h 



-g(t) + — [ dr cos{2u(r - t)}f(r) 
-2uj I dr sin{2w(r - t)}g(r) 

J —oo 



(3.142) 



We further study the correlation among the photons. The time-evolved vacuum state 



\i(t ->■ oo)) = ^2A mfl \m) 

m=0 



(3.143) 



will give 



26 



N = (i(t -> oo)|a i "ajz(t -> oo)) = ^mP m , = |£| 2 - — Jm(/3 2 £), 



m=l 



(3.144) 



AT being the average number of photons in state \i(t — > oo)). Also 



49 



(i(t -+ oo)|aVaa|^ -> oo)) = £m(m- l)P ro , = K| 4 - ^ (M&f) + |£| 2 M&0) . 



m=2 



(3.145) 
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This, to linear order in 9, gives the 2nd order correlation among photons with zero time 
delay to be equal to (see Appendix B.4) 

,»(„) = 1 g /m <« 2 > = 1 - 'iMMl (3 146) 

For the case Im(Ps^ 2 ) < =>- <?^(0) > 1, the photons try to bunch together while for 
Im(/3 3 £ 2 ) > # (2) (0) < 1, they anti-bunch [118]. For the functions (3.79), we get 

r <Rt*\ 2(/ 2 + m^ 2 )sin 2 o;T . 

Im(/3 3 £ ) = ^ 5f sm2wT, (3.147) 

which implies that the bunching or anti-bunching will depend only on the sign of the 
factor g sin2co>T. For the choices taken in figures 3.2, 3.3 & 3.4, uT = | and hence no 
bunching or anti-bunching occurs. 

As a future work one can try to formulate the scattering process in higher dimen- 
sions and study its implications in quantum optics. The correspondence found between 
noncommutativity and quantum optics also encourages one to study such possibilities in 
other forms of time-sapce noncommutativity. As an example one can start with assuming 
the spacetime dependent noncommutative parameter 9 [124-126]. 



Chapter 4 

Thermodynamics of Ideal Gas in 
Doubly Special Relativity 

4.1 Introduction 

In DSR, apart from the constancy of speed-of-light scale, the Planck length l P or equiva- 
lently Planck energy k is also constant under coordinate transformation from one inertial 
frame to another. This leads to modification in the dispersion relation. Consequences of 
the modified dispersion relations on the thermodynamics are being studied extensively 
to infer the effect of Planck scale physics [62-68] . 

The present chapter aims to study the thermodynamics of an ideal gas consisting of 
massive particles in DSR scenario. Both the modification in the dispersion relation of 
the constituent particles and the presence of a maximum energy scale are expected to 
contribute to new effects. DSR transformations can be of several type. We follow the 
formulation of [57, 58] where the modified dispersion relation becomes 

£ i-p2 = m 2 (l-^y . (4.1) 

As < (l — < 1, the energy of a particle with a given momentum decreases in DSR 
with respect to that in SR. This has consequence on the thermodynamics as we will see 
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in § 4.3. The parameter m can be called "invariant mass" as it remains invariant under 
a DSR transformation. Note that in contrary to the SR case, m is no more the rest mass 
energy of the particle. To get the rest mass energy m , we put p = in (4.1). We get 
two expressions for m , namely 

K 

The two solutions are connected by the redefinition of the parameter m — > —m. Hence- 
forth, without any loss of generality we use 

TYl 

™o = — ^. (4.3) 

K 

The physical world is characterized by E < k [58]. In this sub-Planck regime (E p=0 = 
m < k), the positivity of rest mass (m > 0) restricts the range of the invariant mass 
to < m < oo. Thus, in (4.1), we have < p, E < k and < m < oo. 

We study the thermodynamics of an ideal gas in DSR setup. We obtain a series 
solution for the partition function and compute the various thermodynamic quantities. 
We show that our results go to the standard results in the SR limit (k — > oo) [127] as 
well as in the massless DSR limit 



4.2 The Partition Function 

We consider a gaseous system of non-interacting particles obeying Maxwell-Boltzmann 
statistics whose macrostate is denoted by (JV, V, T) where N is the number of particles 
in the system confined in volume V at a temperature T. In the canonical ensemble the 
thermodynamics of this system is derived from its partition function [128] 



Z N (V,T) = J2^M-^E], 

E 



(4.4) 
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where (3 = and ^ denotes sum over all the energy eigenvalues of the system. The 

E 

total energy E of the system can be written in terms of single particle energy e 

E = Y,n £ e, (4.5) 

e 

where n e is the number of particles in the single-particle energy state e and satisfy the 
following condition 

5> e = N. (4.6) 

£ 

We can rewrite Zn as 

Z N (V,T) = ^K}exp[-/3^n £ 5], (4.7) 

where g{n £ } is the statistical weight factor appropriate to the distribution set {n £ }. 
The summation goes over all distribution sets that conform to the above restrictive 
condition (4.6). For Maxwell-Boltzmann statistics, it can be shown [128] 

Zn(V,T) = ^[Z 1 (V,T)] n , (4.8) 

where Z\ is the single particle partition function given by 

Z 1 (V,T) = J2^M-^}- (4-9) 

e 

While for ordinary spacetime, it is easy to show that in the large volume limit one can 
replace the sum by an integral [128] 



Chapter 4. Thermodynamics of Ideal Gas in Doubly Special Relativity 74 



for more exotic spacetimes the measure of integration is expected to get modified 

d 3 p^f(p)d 3 p. (4.11) 

Hence putting together (4.9), (4.10), (4.11) and taking h = k B = 1 we get 

Z 1 (V,T) = -^-3 f d 3 p f(p)ex P i-P(e-m )]. (4.12) 
[2%) J p =o 

Note that in accordance with standard practice, we have subtracted the rest mass m Q 
from the relativistic energy e of the particle. Although there have been few attempts 
[36,67, 129], the form of f(p) is far from settled. Assuming isotropy of spacetime we may 
take f(p) = fip). For a possible deformation of the integration measure, f(p) should be 
expandable in Taylor series in - 



n=0 

with a — 1 since in the limit k — >■ oo, /(p) — > 1. Hence Zi(V,T) becomes 

Z^V.T) = -^3 f rf 3 p E^(^) n exp[-/3(e-m )] (4.14) 

- E^(-o-|) Zf(V,n, (4,5, 

71=0 V ' / 

where (V, T) is the single particle partition function with unmodified measure 

z i (V, T) = -^-3 f d 3 pexp[-/3(e-m )]. (4.16) 

{ZTT) J p =o 

The derivation of (4.15) from (4.14) involves two steps. Firstly, the interchange of the 
summation and integration which is allowed if (see theorem 1.38 of [130]) 

oo i i „ K OO I i / f) \ n 

£^i/^H^-^ = £^(™o-|) zfwrxoo. (4,7, 

TX — Tit — 
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Secondly, writing 



I d 3 p e n exp[-/3 (e - m )] 

Jp=0 



= 4:71 dp p 2 e n exp[-/3 (e - m )] 
Jo 

= 4vr ^m - -^j dp p 2 exp[-/3 (e - m )] 

= ( m °~^) J^d 3 p exp[-P(e-m )] (4.18) 



since the integrand p 2 e n exp[— (3 (e — mo)} remains to be continuous and bounded for 
p £ [0, /c], f3 G [0, oo] (see § 5.12 of [131]). Hence our problem has boiled down to solving 
the integral in (4.16) where e and p are related by the modified DSR dispersion relation 
given in (4.1). The solution of Z\ in the massless case obtained in [69] is 



7 o 2^ 



(2tt) 2 /3 3 



(2 -e-^OSV + 2/3/t + 2)) . 



(4.19) 



The term with e~^ K makes Z\ ml no n- analytic at \ = 0. We anticipate that even when 
m 7^ 0, Z\ continues to be non-analytic at - = and hence does not admit a straight- 
forward Taylor series expansion in K Thus in order to find the leading order deviation of 
DSR thermodynamics from the SR case, one would require a non-trivial series expansion. 



4.2.1 Solution of Z\ 

The dispersion relation (4.1) gives 



p 



2 2 

e — m 



1 - 



K 



1/2 



pdp = 



E+ (1 

K 



de 



Hence changing the variable from p to e in (4.16) we get 



(4.20) 



(4.21) 



Z\ (V, fi) = 



7^2 ex P[/Mi] 



1/2 



exp[— /3e]cfe. 
(4.22) 
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We now consider three different regions of values of m: 



Case I: m — k 



This case is equivalent to the case when the rest mass energy 



K 



m = 



(4.23) 



In this case the dispersion relation (4.1) become 



Also the partition function (4.22) reduces to 



ZUV,T) = 7^ 3/2 expM f K de(2e-K)^e X p[-fie 

[2tt) L 1 J A/2 



(2tt) L ^ J Jk/2 

2y/2V {k\ 3/2 f3 /3k 



(2tt) 2 V/3 



7 



2' 2 7 ' 



(4.24) 



(4.25) 



where 7(0, x) is the Incomplete Gamma Function (see (6.5.2) of [132]). We get a very 
simple analytical form for Z®(V,T). 

Case II: k < m < 00 

This case is equivalent to the case when 



K 



< m < k. 



(4.26) 



We change the variable in (4.22) from e to t — ^ {^Y ~ 1 



— — which gives 



z? 



2\/m 3 



(2vr) 2 [(^) 2 -l 

2^m 3 
(2vr) 2 fe) 2 -l 



3/2 



3/2 



cxp 



cxp 



/3m - 



(3m - 



f3'm 2 

K 
K 



-n/m 



dt t(l-t 2 )^ 2 exp[-f3'mt] 



I*(/3 f m : l)-r (^V-JJ , (4.27) 
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where 



and 



,2 



O -i 



I*(x,y)= I dt t(l-t 2 ) 1/2 exp[-:rt]. 
J-y 

We define Incomplete Modified Bessel function I u (z,y) of order v 



(4.28) 



(4.29) 



such that 



0r"I>+ |) V2 



/i 
(1-i 2 )"^ exp [-zt]dt [Re v>0, \arg z\ < £], (4.30) 



r(x,j/) 



9a; 



7T 



X 



h(x,y) 



(4.31) 



In particular for y — 1, using (3.387 (1)) of [133] and (9.6.26) of [132] we get 



r(z,i) 



ttI 2 {x) 



X 



(4.32) 



where ^(a?) is the 2nd order Modified Bessel function. 

Let us consider a very interesting case of m — > oo (m — > k) 1 . In this limit one 
gets from (4.1) 



-2 J2 



m- 



£- = (!-£ 



J =>. £ — ^ /t Vp G [0, k] 



Thus the total energy E of the system becomes 



E = Nk 



(4.33) 



and the thermodynamics simplifies. Entropy can be computed by counting the total 
number of microstates Qjy available to the system 



N 



N 



2Vk 



3 \ N 



N\ \3(2ttY 



2 I ' 



(4.34) 



1 We would like to thank Diptiman Sen for pointing out this interesting case. 
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where Q 1 is the total number of microstates available for a single particle. Thus the 
entropy S of the system is 



S = In 



2Vk 3 



N\ \3(2ttY 



N' 



(4.35) 



The first law of thermodynamics in this case becomes 



dE = -PdV + ndN. 



(4.36) 



Note that the usual term TdS has been dropped as from (4.35) it is evident that S is a 
function of N and V alone. The pressure of the system is zero as 



P = - 



dE 







N 



while the chemical potential is 



dE 
dN 



K. 



(4.37) 



(4.38) 



v 



Equation (4.16) can now be easily integrated to give 



2V k, 3 
(2^1' 



(4.39) 



which gives the limiting behaviour of I* (/3'm, ^) using (4.29) and (9.6.7) of [132] 



(3m 



K 



(™) 2 _ i m 



(4.40) 



Case III: < m < k 



This case is equivalent to the case when 



< m < 



K 



(4.41) 
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We change the variable in (4.22) from e to t = ^ 1 - (™ ) 2 j + ^ to get 



Z\ = 



2Vm 3 
~ r 9l 3/2 

( 27 r) 2 [l-(f) 2 _ 

2\/m 3 



cxp 



exp 



/3"m 2 



/3m + 



/3m + 



K 



/n/m 
dt t(t 2 - l) 1 ^ ex p [_/3" m t] 



/3"m 



ft 



IT (/Tm, 1) - W /3"m, - 



(4.42) 



where 



and 



(4.43) 



/•oo 

K*(x,y)= dt t(t 2 -l)^ 2 exp[-xt]. (4.44) 
As in Case II, we define Incomplete Modified Bessel function K v {z,y) of order v 

K u {z,y) = r ^i) (|) y exp[-^]dt [ite i/ > \arg z\ < |], 

(4.45) 

(4.46) 



such that 



!T(x,y) 



dx 



Ki{x,y) 



In particular for y — 1, using (9.6.23) and (9.6.26) of [132] we get 



K*(x,l) = 



K 2 (x) 



x 



(4.47) 



where K 2 (x) is the 2nd order Modified Bessel function. 

We shall now obtain the series solution of K*(x,y). We rewrite (4.44) as 



K*(x,y) = J™dt ^(l-^) 



1/2 



-xt 



(4.48) 



Inside the integral t > y and for y > 1 (which is a valid assumption for the case of our 



Chapter 4. Thermodynamics of Ideal Gas in Doubly Special Relativity 80 



1 /9 

interest) the factor (l — ^) can be expanded in series of ^ to get 



with 



and 



where 



K* (x, y) = I dfi t 
'y 



r=l 



dut = t 2 e~ xt dt 



fr(t) 



f 2r' 



(0-|)(l-|)...(r-l-|) _ (2r-2)! 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



r! 2 2r - x r!(r - 1)! 

Now the integral and the summation in (4.49) can be interchanged if (see theorem 1.38 
of [130]) 

OO poo 

J2 dpi t \f r (t)\ < oo. (4.53) 
r =i J y 

Now as t r is —ve for all r > 1 we have 



l/r(0l = -/r(*)- 



(4.54) 



This allows us to interchange the summation and the integral if the final series is con- 
verging. Thus we get 

1 OO 

K* (x, y) = M - -M x + t r M r , (4.55) 



r=2 



if the above is a converging series (see Appendix C.l for convergence of K* (x,y)). Here 



M r 



-f 
J y 



dt t 2 ^e~ xt 



(4.56) 



for r = 2, 3, 
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M and Mi can be easily calculated to be 



exp (—xy) 



x° 



Mi 



((xy) 2 + 2xy + 2) , 



exp (-xy) 



x 



(4.57) 



(4.58) 



Now, changing the variable to t' = xt in M r for r > 2 we get 



oo -t' 
2r-3 I jj e 



M r = x zr ~ 6 / dt 

'xy 



(f/)2r-2- 



(4.59) 



Taking e *' as first function, if we do the integration by parts again and again we finally 
get 



M r 



x 



2r-3 



2r-3 



(2r -3)! 



Ei(xy) + e~*y Yl 



-x 



,fc-i 



^ (2r - 3)(2r - 4)...(2r - 2 - jfe) 



2r-2-fc 



. (4.60) 



Here Ei(x) is the Exponential Integral (see (5.1.1) of [132]). A similar attempt to obtain 
the series solution of I*(x,y) fails! 

Although we obtain the solutions of Z\ in three different regions of values of m, Z\ can 
be shown to be smooth in m (see Appendix C.2) and hence we do not expect any phase 
transition like thermodynamic discontinuity as we vary m. We use the continuity of Z\ 
to obtain the limiting behaviour of I*(/3'm, ^) and K*(/3'm, ^) as m — > k. From (4.25), 
(4.27), (4.32) and (9.7.1) of [132] we obtain the leading order behaviour of I* (0'm, ^) 
as m — ?■ k + to be 



V m/ 



(Bnfl 2 



5 3/2 



e 2 * 



—71 



■ •■ i 3 Bk 
1/2 + 2e"- ' 



7 



2' 2 



(4.61) 



with 8=^-1. For m ->■ re - , using (4.25), (4.42), (4.47) and (9.7.2) of [132] the leading 
order behaviour of K* (/3"m, ^) turns out to be 



K*(P"mA) m ^ 7^372- 



1/2 o A 3 (3k\ 



(4.62) 
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where e = 1 — — 

K 



4.2.2 Leading order deviations 

Having obtained the series solution of Z\ in Case III, we shall now obtain the leading 
order corrections from the massless and the SR cases. 

Leading order deviation from the massless case 

Thermodynamics of a photon gas in DSR with dispersion relation (4.1) and unmodified 
measure has been worked out in [69]. Here we calculate the deviation of single particle 
partition function from that of a photon gas. On expanding Z\ in rj = — with m — > 
(assuming n to be finite) and using (9.6.10) and (9.6.11) of [132], we get 

%\ — Zlml + Zlmlcorn (4.63) 

where Z\ ml is the single particle partition function of photon gas in DSR scenario with 
unmodified measure [69] and Z\ mlcorr is O [r\)\ 

Zlmi = ^(2-e-^V + 2^ + 2)), 
Z^icorr = ~ (2^3 ^ Hv) W + O^)) + (f} K Z° lml ) V + 0( V 2 ). (4.64) 

Note that the correction due to mass of the constituent particle is non-perturbative in 
nature as the first term in Z\ mlc(yrr which contains ln(r^) is the non-analytic piece and 
does not allow a trivial Taylor series expansion at r\ — 0. We can rewrite (4.15) as 

Z\ = Zi„ d + Z lrn i cor ., n (4.65) 

where 

n=0 ' 
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'lmlcorr 



70 

J lmlcorr 



n=l 



n\n' 



8 n 7° 

lmlcorr 



d 



n-l 70 



dp 



•qriK- 



lml 



dp 



n—l 



(4.67) 



The above leading order behaviours have been plotted in Fig 4.1. For our choice of 
parameters they match with the numerical plots up to ^ ~ 0.012. 
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Figure 4.1: The single particle partition function is plotted vs m (left) and - (right). 
The plots for the unmodified measure are in solid lines while those for modified measure 
are in dashed lines. The plots in massless (left) and SR (right) cases are shown in red. 
The leading order behaviours are plotted in blue. The numerical plots for (4.42) and 
(4.15) are shown in green (solid and dashed respectively) for comparison. Different values 
for the parameters in natural units are as follows: V = 10 35 ,iV = 10 25 ,T = 0.01, = 
1, ai = 0.2, a 2 = a 3 = ... = and n = 1 (left), m = 0.01 (right). 



Leading order deviation from the SR case 

On expanding Z\ in rj = ^ with k — > oo (assuming tjiq to be finite), we get 



7° — 7° _L 7° 

^1 — ^ISR ^ISRcorri 



(4.68) 



where Z\ SR is the single particle partition function in SR and Z\ SRcorr is O [r\). 



_ 2Vm| pmo K 2 {pm ) 



'1SR 



(27T) 



f3m 



(4.69) 



70 

J lSRcorr 



I -^H)Prn Z° SRV + O( V >). 



K 2 (f3m ) 



(4.70) 
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Note that the DSR correction is non-perturbative in nature as the first term in Z\ SRcmr 
which contains e -/3m °/' ? is the non-analytic piece and does not allow a Taylor series 
expansion at rj = 0. This is a novel feature in DSR as we know that SR thermodynamics 
is perturbative in the non-relativistic limit: 

±nV(k B Tm f 2 1/2 
Z 1SR - — u e K 2 {u) 





= *i.y W (1 + ^ + 0( (") ) )• (1-71) 



where u = and Zinr is the single particle partition function in the non-relativistic 
case. We can rewrite (4.15) as 

Zl = ZiSR + ZiSRcorr, (4.72) 

where 

Zisr = Z° SR , 
r _ 7 o , Q i m o 7 o fl i ® z isr 

^lSRcorr ~ ^lSRcorr ^ — ^lSR ~ ~ ■ \^-' 6 ) 

The above leading order behaviours have been plotted in Fig 4.1. For our choice of 
parameters they match with the numerical plots up to ^ ~ 0.08. Having obtained the 
leading order correction to Z\ due to DSR, we shall now compute its effect on the various 
thermodynamic quantities. 



Chapter 4. Thermodynamics of Ideal Gas in Doubly Special Relativity 85 



4.3 Thermodynamic quantities 

The free energy F, pressure P, entropy S, internal energy U, internal energy density p 
and heat capacity Cy are defined as 

F = -^lnMAmo)) = -^|ln^j+l} (4.74) 



U = F + TS (4.77) 
P = \ (4-78) 

c * - (i),. 

The above quantities for Case I can be found by using (4.25) and they have been plotted 
in Fig 4.2. Now we shall obtain the leading order thermodynamics for Case III. If we 
denote the free energy, pressure, entropy, internal energy, internal energy density and heat 
capacity obtained in the SR or massless cases by F , P , S , U , p an d Cyo respectively, 
and write Z\ = Z w + Zi corr , where (Z w , Zi corr ) = (Zisr, Z-^sRcorr) or (Zi m i, Z\ m i corr \ we 
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+ 
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(4.80) 



(4.81) 
(4.82) 



lcorr 



z 2 w dp dp 



no 



d 2 P 



n dZi co , 
Zio dp 



+ o 



Z\ t 



no 



(4.83) 
(4.84) 



where n = y is the number density. The correction to F depends on the ratio of Zi corr 
and Zio and is independent of the volume V of the system. Note that this is true to all 
orders. Hence, the pressure P of the system which is defined as P = — (^) NT gets no 
correction: 

P = Po. (4.85) 
The equation of state in SR is (see (8.128) and (8.134) of [127]) 



SR 



PSR 



which gives the following DSR equation of state 



P + 



~~P + 



n dZisR 

C 

Zisr 9/3 



im + 



K 2 {0m) 



Pm 



+ o 



ZlSRcorr \ 
ZlSR J 



(4.86) 



(4.87) 
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In Fig 4.2, we have plotted the various thermodynamic quantities for Cases I and III 
as a function of T and compared them with the SR case. The P vs p plots have been 
obtained by varying T keeping all other parameters fixed. The qualitative natures of the 
plots for different cases are same. In case of S there are two competing effects: while 
the cutoff tries to reduce S by limiting the number of accessible states, the modified 
dispersion tries to increase S by enhancing the Boltzmann weight exp(— e/T) (note that 
£dsr(p) < £sr(p) for a given momentum state p, the change being more for greater value 
of the parameter m). At low temperatures, the latter is dominant and Sdsr > $sr- For 
our choice of parameters this is clearly visible for the plot of S in Case I. In the high T 
regime, the cutoff effect comes into play and Sdsr < Ssr- The cutoff also saturates U 
as T increases, and Cy — > 0, resulting in a steeper equation of state. Here we make an 
interesting observation. There have been attempts to define velocity in DSR [134]. If 
we adopt the usual definition for the speed of sound c s = then we observe that c s 

grows without any bound. Possibility of such scenarios has been discussed in [135]. 

For given choice of parameters in case of quantities like S and U, the leading order 
behaviours for Case III match with the numerical plots up to T ~ 0.15, while in case of 
Cy which contains second order derivatives of the partition function with respect to T, 
the leading order behaviours match with the numerical plots up to T ~ 0.09. Note that 
the leading order behaviours have been obtained assuming to be finite and — — > 
which in turn implies ^ = ^ — > 0. Hence as T increases, the leading order plots depart 
from their numerical counterparts. 
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Figure 4.2: The plots for the unmodified measure are in solid lines while those for 
modified measure are in dashed lines. The SR plots are shown in red. The figures in the 
left column are for Case I while those in the right column are for Case III. The DSR plots 
for Case I are plotted in green. For Case III, the leading order behaviours are plotted in 
blue. The numerical plots for Case III are also shown in green for comparison. Different 
values for the parameters in natural units are as follows: V = 10 35 ,iV = 10 25 , k = 
1, m = 1 (Case I), 0.01 (Case III), a — 1, a x = 0.2, a 2 = a 3 = ... = 0. 



Chapter 5 



Conclusions 



In chapter 2 we described static classical solutions of noncommutative gauge theories in 
various spacetime dimensions and showed that the GBOs are significant in constructing 
solutions with higher topological numbers. We started with describing the static mag- 
netic flux tube solutions of higher moment found by Polychronakos in terms of GBOs. 
In doing so we get a better understanding of the winding numbers n j's. They corre- 
spond to the radii of the flux tube in different irreducible parts of the full Fock space. 
The increase in the degrees of freedom (higher number of required noi's) due to use of 
reducible representation of the oscillator algebra cause to expand the solution space of 
the flux tube solitons. On the other hand, the vortices with higher winding numbers 
correspond to known solutions. 

The case of multi-instantons is different. The multi-instantons with charge —p\p 2 
(pi,P2 non-negative integers) are not gauge equivalent to known solutions. Another 
significant result of this study is an explicit relation between the instanton number and 
the representation theory labels pi and p 2 - The charge of a multi-instanton does have 
the information about the reducibility of the space along with the topological nature 
of the solutions. Using the "translated" b operators (2.41) we have constructed multi- 
instantons that depend explicitly on p\p 2 complex parameters. While the full moduli 
space of noncommutative multi-instantons is still not well understood, we hope that this 
identification contributes partially to this question. 
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Though we have only considered a few cases, there is actually a large variety of 
situations in noncommutative gauge theories where GBOs may be used. In particular we 
expect that this procedure may shade new light on merons, monopoles, dyons, skyrmions 
etc. We plan to revisit some of these questions in future. 

In chapter 3 we studied the effect of noncommutative time coordinate. To do so we 
took a simple quantum mechanical system with time-space noncommutativity in 1+1 
dimension. We developed a formalism to compute the transitions between asymptotic 
states of the quantum mechanical system with noncommutative time. We found that 
for a free Hamiltonian in R^' 1 which is independent of time, the transitions are equal to 
the same for a different Hamiltonian in R 1,1 found after the replacements (4.10). The 
time evolution of an operator and its expectation value (and hence also its uncertainty) 
can also be found in a similar manner. Specifically, for FHO the transition probabilities 
get modified and is given by (3.64) and (3.72). The Poissonian distribution for the 
"vacuum to any state transition" also gets modified and is given by (3.77). The study 
of uncertainties in position and momentum says that the time-evolved state is no more 
coherent. It gets some squeezing effect due to the noncommutativity, keeping the product 
of the uncertainties minimum. These uncertainties are explicitly found and is given in 
(3.127), (3.128) and (3.129). The leading order corrections in these uncertainties are 
oscillatory in time and they depend independently on the mass of the particle m and 
the frequency of the oscillator u (note that the commutative uncertainties depend only on 
the product m u). These results suggest a possible modelling of the noncommutativity 
for the nonlinear phenomena in Quantum Optics. The noncommutativity results in the 
following nonlinear effects: 

1. The photon-count gets modified from the usual Poisson distribution. 

2. The uncertainties in the field quadratures change keeping the product minimum 
(the squeezing effect). 

3. The second order correlation function g^ 2 \0) gets modified producing new effects 
like bunching or anti-bunching of photons depending on the value of im(/3 3 £ 2 ). 
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All these observations suggest that the noncommutativity produces incoherency in the 
otherwise coherent field. 

In chapter 4 we used the dispersion relation (4.1) and have considered the modified 
phase space measure (the modification being isotropic and expandable in Taylor series). 
We considered three cases separately (m = k, m > K,m < k). The single particle 
partition function has been shown to be smooth in m G (0, k) (see Appendix C.2). For 
the case m — k, a simple analytical form for the partition function has been obtained (see 
(4.25)) while a series solution for the partition function has been obtained for m < k (see 
(4.42) and (4.55)). In doing so, new type of special functions (Incomplete Modified Bessel 
functions) emerged. We observed that DSR thermodynamics is non-perturbative in the 
SR and massless limits. Using the leading order solutions, we derived thermodynamic 
quantities like the free energy, pressure, entropy, internal energy and heat capacity (see 
Fig.4.2). A stiffer equation of state has been found. 
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A.l Computing the flux tube solutions 

We start with the ansatz 

V = af(N), V = -V ] = -f*(N)<f. (A.l) 
The equation of motion (2.66) reduces to 

af(N) {(N + l)\f(N + 1)| 2 - 2N\f(N)\ 2 + (N - l)\f(N - 1)| 2 } = 0. (A.2) 
We write f(N) in the basis of number states 

oo 

f(N) = J2f(n)\n)(n\ (A.3) 

n=0 

to get 

/(n) = 0, or (n+l)|/(n+l)| 2 -2n|/(n)| 2 + (n-l)|/(n-l)| 2 = 0. (A.4) 

We try to find a solution for which f(n) satisfies the first of the above equations, i.e., 

f(n) = for < n < no and the second of the above equations for n > n . This choice 
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leads to the localized flux tube solutions as discussed in the section 2.3.1. The most 
general solution of this type is given by 

for < n < n 

fin) = { / " " (A.5) 

■ e^"|/(n + l)h/ (no+1)(n ~ no) forn>n 



for < n < n 

\f(no + 1)| V / (no+i ; n -" o) for n > n 



For the choice (f n = 0, \f(n + 1)| = , 1 we get 

V (n-o+l)0 



n=no 



(A.6) 



which corresponds to the solution (2.68). 
Let us now start with the ansatz 

V = bG{N). (A.8) 

The equation of motion reduces to 

bG(N) {{M + 1) \G(N + 2)| 2 - 2M |G(A^)| 2 + (M - 1) |G(A - 2)| 2 } = 0. (A.9) 
We expand G(N) in the number state basis, i.e., 

oo 

G(N)=^2G{n)\n){n\ (A.10) 

n=0 

to get 

G(n) = 0, or Q + l) |G(n + 2)| 2 - n|G(n)| 2 + Q - l) |G(n - 2)| 2 = 0. (All) 
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for n = even and 

G(n) = 0, or + l) |G(n+2)| 2 -(n-l)|G(n)| 2 + - l) |G(n-2)| 2 = 0. 

(A.12) 

for n = odd. There exists a solution of the form 

for < n < nn + 

G(n) = { , " " (A.13) 

' |e^+" \G(n 0+ + 2)| x / ("o + +2Xn-n 0+ ) f or n > nQ+ 

(A.14) 



2 

for < n < n 0+ 



e*+» |G(n 0+ + 2) | ,/ (no++2)( "~" ()+) for n > n 



for n = even and 



for < n < nn 

G ( n ) = <,... ., /t- : ; " " (A.15) 



ie^-» |G(rao- + 2)| y / (no - + ^-"°- ) for n > n _ 
for < n < no- 



§e*-» \G(n _ + 2) | / ("°-+iX"-"°-) for n > ^ 



(A.16) 



for n = odd. Here n + and n _ are even and odd integers respectively. For the choice 
G(n 0+ + 2) = 2 G(n - + 2) = 2 y? +n = ¥>-„ = we get 



G(iY) = -L^/ JV -"° + ^-"°- A - ^ |n>(n|A ++ £ |„) <n |A_j (A.17) 
which corresponds to the solution (2.73). 



oo 



A. 2 Determination Of and 

We multiply to the first equation of (2.93) from right and $oo to the second equation 
of (2.93) from left and use (2.90) to get 

4 = -i(a- $ooa$oo) , A 00 = i(a j < - Q^Q^) . (A.18) 
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For the solutions (2.91) we get 



= -i (a - ^ood^oo) (A.19) 

= -|- -^ + ")' Af + "- 1 ' <" + 1 >'| (A.20) 
y/(N + n-l)(N + n-2) N 1 



-ia j= (A.21) 

y/N 



= ~ i 7WTt a ^~^ wT ^y (a,22) 

Here we have used the identity 

aV n = (N + n)(N + n-l) (N + l). (A.23) 

In a similar fashion we get the expression for A^. The same technique can be used to 
compute A™ w and A™ w as in (2.99). 



A. 3 Large Distance Behavior 

Let \u) be the coherent states of the operator a, where a; is a complex number labeling 
the states. In the large oj limit, the expectation value (oj^^oj) and (^A^u) gives the 
large distance behavior of the solution: 



= (cu\ a n \co) (A.24) 
Va n a' n 

= cu n (cu\ , 1 \u) (A.25) 

« u n (u\ j_ \u) (A.26) 



« B (H^I«» 5 (A.27) 
w n -^r (A.28) 
e ir T (A.29) 
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For the gauge field, 



(VN - VN + n)\co} (A.30) 

» -i(u\ a (l-(l+^L)|q;) (A.31) 

~ igMa-la;) (A.32) 
ra 1 

<£. (A.34) 



Similar thing can done for and A™ w to get (2.101). 



A. 4 Computation of the leading order magnetic field 
for the Nielsen-Olesen Vortices 

Using (A. 18) in the expression of (2.88) we get 

Boo = -1+ [$00^00, $00^00] (A.35) 

= $oo (al» 0O $ 0O a t - a^oofcood) $00 - 1 (A.36) 

Using (2.92) one can calculate 

n-2 

a^oo^ooa 1 " - a^oo^ooO = 1 + n|n)(n| + (ra - l)|n - l)(n - 1| - ^ |m)(m|. (A.37) 

m=0 

Finally for the solution (2.91) we get the magnetic field as given in (2.97). Similar 
technique can be used to calculate B^ w . 
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A. 5 The New Solution Satisfies The Equation Of 
Motion 

The equations of motion are 

[D xP ,[D xP ,D x a]]=0, a, (3 = 1,2,3,4. (A.38) 

There are four equations for different values of a. We will show only for a — 1 and other 
follow similarly. For a = 1 the equation becomes 

[D x2 , [D x 2,D xl ]] + [73,3, [73,3,73,1]] + [73,4, [73,4,73,!]] = 0. (A.39) 

In terms of the complex coordinates, 

([73!, [73!, 73 x ]] + [73 2 , [D- 2 , 73 x ]]) + ([D- u [D- u D,}} - [D- 2 , [73 2 , D-^ = O. (A.40) 

Now using the explicit expression for the solutions (2.164) we get 



, t /A7 + 3 / A7 f M • 2 M 1 , , M M + 3 



[CICCJ] = S^/— ^— I - — M, - — — (M, + 1) 



M M + 3 
M + 2M+ 1 



(M ' + 1) + ^rl^( Ml + 2 )) s - (a ' 41) 



, A rA A .. „ I+ /A7 + 3 / M /M + 2M-1„ M M + 3„ 

ID,AD,D,]] = -S.„ Jl ^ — ^ — ( — — — - M 2 - — — - t M 2 

A7 M + 3 
~A7 + 2 A7 + 1 

Adding above two we get 



01, [£>l,Di]} + [D 2 ,[D 2 ,Di]] = 0- (A.43) 
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Similarly we can show [Di, [_Dj, Di]] — [D 2 , [D2, Di]] = and hence the equation of motion 
is satisfied. 

A.6 The New Solution Satisfies The ASD Condition 



The commutator between D a and D a is 



1 



[D a ,D s \ = --S, 



M M + 3 



9 KM + 2M + 1 



(M a + I) 



M + 2M- 1 



M M + 



-M a \Sl w . (A.44) 



Summing over a = 1, 2 we get 



[Di,Di] + [D % M = -\s new (\r][ • \ - (M + 2)^-p) Sl w = -\. (A.45) 



e 



M + l 



M + l 



e 



Therefore, F ± i + F 22 — § + [-Di, -Di] + [A2, -D2] — 0- Similarly one can show F 12 = F\ 2 = 0. 



A. 7 The Topological Charge Of Instantons 



The topological charge of the instanton is defined as 



Q^9 2 Tr (F n F 22 - F l2 F 2l ) Q = c9 2 ^ (m, n 2 \ {F n F 22 - F l2 F 2l ) |m, n 2 ) 

ni,n2=0 

(A.46) 

where we denote the constant of proportionality by c. 
Further, the ASD condition gives 



/ 'lJ / '22 /fl - ~ ( 02 + g 



+ 



(A.47) 
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A. 7.1 Usual Single ASD solution 

The solutions (2.146) can be simplified to 



£>*=±sJ f /( N + 3) s a a S\ * = -±Sad"( N + 3 \,St. (A.48) 
V(N + l)(N + 2) a ' ^6 a \{N + l){N + 2) 1 ; 



An arbitrary commutator becomes 

Da,D~ b 



-ScJa 1 gt _ *<* s + t fA m 



and the components of the field strength (and their products) can be calculated to be 

f-f-- 4 9 {N2 - Nl)2 ?t r A ^m 

FllF22 - ~^V(iV + l)2(iV + 2)^ • (A - 50) 
Again using (A. 49) we get 

F * = -t S<iai N(N+l)(N + 2/> F2l = "e Sa{a2 N(N+ 1 l)(N + 2) Sl (A ' 51) 
Multiplying the above two and simplifying gives 

F - F -- 16 o (N 1 + 1)N 2 

hl2 * 21 ~ ¥ b nW+WW+W ' { ' 

Thus we get 

= _4cVV (n 1 + n 2 ) +4n 2 ^ 1 

~ ^on^i(^ + ^) 2 K + ^ + mn 1 + n 2 + 2y 2^ o{ni + 2 y(n l + 3y 

(A.53) 



Appendix A. Appendices for Chapter 2 



100 



A.7.2 The New Solution 

Again the solutions (2.164) can be written as 



~ 1 / M(M + 3) f h _ 

Ua ~ ^Q bnew Y(M+l)(M + 2) a new a ~ ^9 



Then the commutators evaluate to 



l -s bul M(M + 3) - gt 

new a V( M + 1 )( M + 2 ) 



(A.54) 



= _4 f 1 _<U M(M + 3) t 

6 a M(M + l)(M + 2) ne ™ neu, (M + l)(M + 2) 



(A.55) 



The product of the field strength becomes 



F p -_!<? (M 2 - MQ 2 f 

11 22 ~ M 2 (M + l) 2 (M + 2) 2 neu) 1 ' 



F F - 16 9 ( M i + at 

Hence the charge is 



L2-2J = 02^» M 2( M+1 )2( M + 2 )2' 



(A.57) 



oo / oo 



Qnew 4c ^ ^ j ^ ^ 



(m ni +m„ 2 ) 2 + 4m n2 



ni=0 \ri2=2 



,m ni + m n2 ) 2 (m ni + m„ 2 + l) 2 (m ni + m„ 2 + 2) 



2 



(m ni +2) 2 (m ni +3) 2 

with 



(A.58) 



1 _ _ J ;n = even 

™n = 7}( n -K)> A n = < ( A - 59 ) 

2 \ 1 ;n = odd 

Now for all even n's we have m n = m n+1 . Hence any absolutely convergent series 
over n\ (or n 2 ) whose terms depend only on m ni and m n , 2 can be broken into equal sums 
to give 

oo 

^G(m m ,m n2 ) = 2 G(m ni ,m n2 ). (A.60) 

m=0 ni=0,2,4,... 
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Hence 



oo oo 



ni=0n 2 =2 m=0,2,4,- "2=2,4,6,... 

= 4 ^ F(m ni ,m n2 ). 

ni=0,2,4,...n 2 =2,4,6,... 



(A.61) 



Again m n = | for all even n's. Thus the charge becomes 



Qr 



ni=0,2,4,...n 2 =2,4,6,... I 2 2 ^ I 2 2 1 J I 2 + 2 + Z J 



E E 

=0,2,4 

-4c £ 



(f + f) 2 + 4(f) 



"1=0,2,4,. 



(f + 2) 2 (f + 3)' 



(A.62) 



Redefining ^ — >■ ni and ^ — >• n 2 , we get 



Qnew 4:Q- 



(A.63) 
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B.l Orthonormality of the eigenfunctions of the non- 
commutative SHO 

The eigenfunctions of the noncommutative SHO Hamiltonian Hq given by (3.14) are 
given in (3.16). If we define <f> n (ki) such that 



(B.l) 



we can write 



ij n (x ,x 1 )= / dk dki<p n (ki)5 



kn + ou { n + 



(B.2) 



Hence 



^ n5 (x ,xi) = J dk dk x 4> n {ki)5 



k + to ( n + 



gifcirrigifcozo 



(B.3) 
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4> n (xi) being the orthonormal eigenfuntions in the commutative case 9 = 0. The inner 
product of two eigenfunctions become 



(^ n (x ,x l ),^ m (x ,x 1 ) t = J dx 1 il}* s (t,x 1 )il} m s(t,x 1 ) 

= J dx^lix^^^rnix^e-^^y 

= e'^-^J dx l <Pl{x l )<p m (x 1 ) 



with 



Thus 



(B.4) 



B.2 The Green's operator function 

Expanding G(t,t ) as a Fourier integral we get 

i r+oo 

G(t, t ) = G{t - 1 ) = — dJgiJ)e-^^\ (B.5) 

27r y_ 00 



g(u') = (W - /Jo)" 1 = £ fc , , e-M (U (B.6) 

00 / 1 r+oo 1 \ 

G(t,t ) = Mim - / e -^) (B.7) 

Here, ie has been introduced to avoid the pole on the contour (real axis). After finding 
the integral inside the summation using the complex analysis we get the Green's function 
to be 

G{t,to) = - l T Q{t-to)e-^- t »\ (B.8) 
where Q(x) is the Heaviside Step Function. 



Appendix B. Appendices for Chapter 3 



104 



B.3 Simplifying A^(T 2: T x ) 

To simplify (3.45) what we do is to find the differential equation for A^\t, Ti) with inde- 
pendent variable t and solve it with proper initial conditions. The differential equation 
has been found to be 

(ih^ - H%(t)j AW(t,rO = Hfl^A^^T,), (B.9) 

which is a first order equation and hence A^(t,Ti) is unique if an initial condition 
is given. The initial condition comes from the fact that A^°\t, Ti) and AW(i,Ti) must 
become identity and zero respectively for t — Ti, i.e., no interaction. Thus, A^\Ti, Ti) = 
0. To solve the equation we define the Green's operator function Gi nt (t,t ) as 

ih%- - Hjf(t)) G int (t,t ) = 5(t - t ). (B.10) 



dt 

Now generalizing solution (B.8) for the time-dependent case we get 



G ftrt (Mo) = -!e(f-fo)T e-ww =~e(t-t )A^(t 1 t )- (B.11) 



ft 



h 



It can be easily checked that the above expression for the G int (t,to) satisfies the corre- 
sponding differential equation. The solution for A^\t,Ti) is then given by 

/+oo 
dt GUMo^^oM^o,^), (B.12) 
-oo 

where A^ m (t, Ti) is the solution of the homogeneous equation 

<^- J H)S t (0)^£L(*,r 1 ) = 0. (B.13) 

The O-function in the expression of Gi n t(t,to) and the fact that interaction was off 
before Ti, with the initial condition for A^(Ti,Ti) = 0, gives the initial condition for 
A^ m (t, Ti), i.e., ^^(Ti, Ti) = 0. The only solution of the homogeneous equation with 
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this initial condition is A^ m (t, Ti) = 0. Thus we get 

/+oo 
dt G int (t,t )H}f (^^(fcTx). (B.14) 
-oo 



Now, putting the expression of G int (t, t ) above and introducing A^°\to, Ti) [yl( )(i , Ti)] 
before i?}"*(i ), we get (3.46). Here we have also used the following property of A^°\t, Ti): 

^ft.tojA^ticT!) = A(°)(*,Ti). (B.15) 

B.4 The correlation function 

An operator corresponding to the detection of a photon by a detector should be propor- 
tional to the annihilation operator a (say ka) [118,136]. Hence if \i) is the initial state 
of the radiation field, the state after the detection of one photon is ka\i). The amplitude 
for going to the final state |/) is given by k(f\a\i). The corresponding probability is 
|/c| 2 |(/|a|i)| 2 . Thus the probability of detection of one photon in the state \i) 

Pi = J>| 2 |(/|a|i>| 2 = \k\ 2 J2(i\^\f)(f\*\i) = \k\ 2 (i\^a\i). (B.16) 
/ / 

Similarly, probability of detection of two photons with a time delay of r is 

P2 = 5>I 4 K/I°(* + r)am)\ 2 = \k\*J2(^(tW(t + r)\f)(f\a(t + r)a(t)\i) 
f f 

= \k\ A {i\a){t)a)(t + r)a(t + r)a(t)\i). 

(B.17) 

The 2nd order correlation function with a time delay r is defined as 
(2) P 2 (r) (i\ri(t)a}(t + T)a(t + T)a(t)\i) 

9 (r) = ~W = pj#F • (R18) 

For r = 

(2) (i\d(t)d(t)a(t)a(t)\i) (i(t)\a)(Jaa\i(t)) 
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For a coherent state it can be calculated to be equal to 1. 



Appendix C 

Appendices for Chapter 4 



C.l Convergence of K* (x, y) 



K* (x,y) given by (4.55,4.57,4.58,4.60) is convergent if the following two series are con- 
vergent: 

~ (_ T \2r-3 

S = E'.7#W (CD 



and 



2 (2r-3)! 



oo 2r-3 ( _ x \k-l /^\2r-2-k 



S ' 5' r S(2r-3)(2r-4)...(2r-2-fc)y " (C ' 2) 

5 can be easily proved to be absolutely convergent using Ratio test. For S' first consider 
the following double series: 

oo oo , \k-l /l N2r-2-fe oo oo 

S " = y tr y- — — { ~ X A { - — — -(-) =YYa rk . (as) 

Let us first test the convergence of S" (See theorem (2.7) of [137]). The row series S r 
(for a fixed r) and the column series Sk (for a fixed fc) are defined as 



S r = J2 a r,k, (C.4) 



fc=i 
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Sk — ^ ] a r,k- 



(C.5) 



r=2 



The ratio tests for S r and S k show that they are absolutely convergent (for y > 1). Also 



lim 

r,k— »oo 



< 1. Hence, 5"' is absolutely convergent. Now 



oo 2r— 3 



\ S "\ = '}2^2\ a r,k\ = '}2^2\ a r,k\+ L ]L>0. 



(C.6) 



r=2 fc=l 



r=2 fc=l 



As IS"' | is convergent (or in other words S" is absolutely convergent) we must have 

oo 2r— 3 

\S'\ — ^ y^l^r.fcl to be convergent (or in other words S' to be absolutely convergent). 

r=2 k=l 

Thus the series expansion of K* (x,y) is absolutely convergent. 



C.2 Continuity and differentiability of the partition 
function in mo 

We shall show that is continuous in m for m G [0, k]. After integrating over the 
angular coordinates (4.16) gives 



Z\ (m ) 



2V 
^fJo 



dp p e 



2 p -/8[e(p,mo)-mo] 



2V [ K 
—2 / dp f{p,m ). 
(<Zn) Jo 



(C.7) 



The integrand f(p,m ) = p 2 e -l 3 Hp^ m o)-mo} j g a continuous bounded function of p and 
m in the range m G [0, G [0, k\. Thus Z° (m ) is a continuous function of m 
as the function g{p) = C, where C is the upper bound of \f(p,mo)\ in the range mo G 
[0, k],p G [0,k], satisfies g{p) > \f(p,m,o)\ for all m G [0,«],p G [0, /c] and is integrable 

r 

as / dp g(p) = Ck < oo (see Lemma 1 in § 5.12 of [131]). 
Jo 



The derivative of the integrand with respect to m is given by 



df{p,m ) 
dm 



= Pf(p,mo) 



\ 



e 2 (p, m ) 



pf 



«o (1-^) 1/2 (p 2 



-i 



1/2 



(C.8) 
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It has 2 poles (and also branch points) in the complex p-plane at p — ± ^ 2m ™° y /2 . We 
note that the poles and the branch points remain to be at the same positions for all 
higher order derivatives of f(p,m ) with respect to tjiq. For m = both the poles are 
at p = and as m Q increases the poles separate towards the imaginary axis. They keep 
on moving on the imaginary axis till they reach ±ioo at m = |. After that they start 
to come closer to each other on the real line and finally at m = k they stop at p = ±k. 
Note that for all m e (0, k) the poles are never on the contour of integration (the 
real line from p — to p — k) and the functions ( ^^™\ 9 g^r^ ? etc.^ remain to be 
bounded. This (by the same argument as given in the case of the continuity of Z\ (mo)) 
ensures the infinite-order differentiability of Z\ (m ) in m G (0, k) and the derivatives 
can be found by using the Leibniz rule (see Lemma 2 in § 5.12 of [131]). Note that the 
fact that we are not being able to say about the differentiability of at m = could 
be a relic of the non- analytic part in (4.64). 
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